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Interferometer for Large Surfaces 


James B. Saunders and Franz L. Gross! 


(n interferometer is described that permits the testing of large areas, such as layout 


plates. 


reflect from the specimen at a large angle of incidence. 


The extension to large areas is obtained by causing a collimated beam of light to 


The resultant fringe pattern is a 


contour map of the surface relative to an arbitrarily chosen plane and the contour interval 
is a function of wavelength and angle of incidence. 


1. Introduction 


An interferometer for testing large surfaces has 
been described by Linnik.? However, in spite of 
certain advantages, the Linnik interferometer has 
not been used extensively in this country. This 
lack of application may be due to the complexity 
of the instrument. 

The instrument described here is relatively simple 
and easy to operate, is twice as sensitive as the 
Linnik interferometer, and is relatively free from 
vibrations. 


2. Optics of the Interferometer 


The optics of this interferometer are shown 
schematically in figure 1. The light from a source 
at S is collimated by the lens 1;, and separated into 
two coherent beams by the semireflecting dividing 
plane of a Kdésters double image prism.** The 
entrance angle of the collimated light may be 
adjusted to give any desired deviation @ between 
the component beams 1 and 2. Beam 1 is reflected 
normally from mirror .\/, and returns upon itself. 
Beam 2 is reflected from the surface to be tested at 
an angle of incidence of 90—8 degrees, then nor- 
mally from mirror .\/, and returns along its previous 
path to the dividing plane of the prism, where it 


a 


P vidress: University College, London, England 
V. P. Linnik, Compt. rend. acad. sci R.S.8. 35, 16 (1942 
J. B. Saunders, Construction of a Késters double-image prism, J 
N Bs, 58, 21 (1957) RP2729 
‘The Kdésters double 
of which is partially silvered ¢ 


ire then cemented together to for 
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prism is made from two 30°-0°-90° prisms, one 
m the face opposite the 60° angl Che two prisms 
m the equilateral prism shown in figure 1 Phe 
part v silvered surface becomes a semireflecting plane 


recombines with beam 1. The observer at £ sees 
interference fringes on the superimposed images of 
the two mirrors. 

The gross aspects of the fringe pattern (fringe 
direction and spacing) are controlled by the wedge 
angle between the wave fronts of beams 1 and 2, while 
the small irregularities of the fringe pattern are a 
function of the irregularities of the test surface. If 
the test surface were perfectly flat, the fringes would 
be straight and parallel. Any curvature of the test 
surface introduces a corresponding curvature in 
the wave front of beam 2, and this wave front, 
when compared with the plane wave front of beam 1, 
introduces curvature into the otherwise straight 
fringes. By adjusting the tilt of the test surface 
(i.e., setting B=@ and introducing a slight tilt 
across the width of the test area), it is possible to 
adjust the fringes so that they run parallel to the 
long dimension of the test area (as seen by the 
observer at /). In this case, the curvature of the 
fringes is a direct and precise measure of the curva- 
ture of the test area (in its long dimension). The 
fringe pattern is a contour map of the area of the 
specimen surface that is being tested, 


3. Sensitivity 


The sensitivity of the instrument depends upon 
the value of 6. Since the light is reflected twice 
from the specimen surface, the sensitivity is double 
that obtained. with the Linnik interferometer for the 
same angle of incidence. When the instrument is 
adjusted as described in the preceding section, one 
fringe departure from straightness corresponds to a 
departure from flatness of \/(4 sin 8), where X is the 
wavelength of light. 


SURFACE TO BE TESTED 


FIGURE | Optica a 


randement of the surface 
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j 
plate interferometer 








In a finished instrument, the angles @ and 8 could 
be made adjustable so that any length surface could 
be made to fill the aperture of the system. This 
would give maximum sensitivity for any size surface 
measured However, since the value of 8 must be 
known for evaluating the fringe pattern, it may be 
more practical to use fixed values for @ and 8 (with 
M, and M. bound into a rigid unit) and bind the 
mirror unit to the prism housing. This would elimi- 
nate the necessity of making frequent measurements 
of 8 and also add stability to the instrument. The 
resulting instrument would have a fixed adjustment 
and could be used by an unskilled operator. 

The maximum length of surface that may be cov- 
ered with one setting is .1 ese 8B, where A is the wper- 
ture of the prism. Thus by decreasing 8, any length 
surface could be covered with a prism of a given 
aperture, but of course the sensitivity of the instru- 
ment would also decrease. 


4. Experimental Model 
Figure 2 is a photograph of the pilot model that 


was used to test the interferometer. The surface 
plate used in this assembly was 91 em long. The 








Pilot model of surface plate interferometer for testing 


FIGURE 2 
optical principals with a typical granite 


feasibility of the 
surface plate 
to give pproximately equal 


designed 


The reflectivity of the end mirrors is 
intensity for the two lightbeams at the receiving point E 


wedge angle between the wave fronts was controlled 
by adjusting the end mirrors, which produced the 
same effect as tilting the surface plate. The Késters 
double image prism and lens assembly shown in 
figure 2 were taken from another instrument that 
was designed for use with a clear circular aperture of 
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51-mm diam (2 in.).° Since the prism had an aperture 
5l-mm square, better coverage of the specimen sur- 
face could have been obtained by increasing the 
circular aperture to 5ly2 mm. This would have left 
the square aperture of the prism unobstructed, so 
that the area of the specimen being tested would 
have been a rectangle 51 mm by 91 em, and the un- 
favorable narrowing of the field at the ends of the 
specimen due to the circular aperture would have | 
been eliminated. 


5. Results | 


The interference fringes shown in figure 3 were | 
taken with the instrument described above and | 
shown in figure 2. Photograph A was made with 
white light and the others with monochromatic 
\=5876 A) light. The difference between photo- 
graphs C and D was obtained by changing the angle 
between the two interfering wavefronts, thus chang- 
ing the width of the fringes. Photographs A, C, and 
D of figure 3 were made with the plate of black 
granite, shown in figure 2, whereas photograph B 
was made with an old cast iron plate (dated 1918) 
that was badly scratched and marred. In this model 
the angle 8 was 2%° and thus a departure from 
straightness of one fringe corresponds to a deviation 
from flatness of 5.75 X (approximately 0.00013 in.). 

The — Hees lod, e of this instrument was 
found to be highly satisfactory. Some mechanical 
deve tht nt is necessary in order to produce a more 
practical instrument. Accordingly, a few suggestions 
given here may be of help to the designer. The prism 
should be designed to use its entire rectangular aper- 
ture. The apertures of the lenses should be suffi- 
ciently large to prevent constriction of this aperture. 
The prism and end mirrors should be rigidly bound 
together if vibrations of the fringes are to be avoided. 
White light may be used but the resulting difficulty 
of finding and adjusting the fringes exceeds the in- 
conveniences associated with monochromatic sources. 





ned for use with @—0 Using it at #2! 
cular aperture became an oval aperture with a vertical diameter of 45 mm, | 
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FIGURE 4. 


Modified Koésters double-i mage prism for use with 
monochromatic light. 


A cover that encloses the light beams, between the 
prism and end mirrors, greatly enhances the stability 
of the fringes, and should be used unless the surround- 
ing air is in a steady state. If white light is to be 
used, differential refraction may be reduced to a 
minimum by using the form of Késters prism shown 
in figure 4. The procedure for making such a prism 
differs from that described in footnote 4 in that all 
surfaces are finished before cutting the prism into 
two parts. 

December 8, 
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Glass Formation in Polymers: II. The System 
Rubber-Sulfur 


Gordon M. Martin and Leo Mandelkern 


Volume 


composition 


versus temperature 
have 


relations encompassing a 
been determined for a series of 


wide range in temperature and 


natural rubber-sulfur vuleanizates and 


mixtures. The characteristic temperature of glass formation is deduced from these studies. 
The glass temperature is found to increase continuously from O° to 90° C in the vul- 
canizates as the amount of bound sulfur is inereased and to be invariant with the composi- 


mixtures 
theory of glass formation. 


tion in the 
free-volume 


1. Introduction 


The physical and mechanical properties of bulk 
amorphous polymers depend to a large extent on the 
relationship between the temperature of study and 
the temperature at which formation 
Thus at temperatures above that for glass formation 
amorphous polymers exhibit rubberlike elastic be- 


glass occurs. 


havior, while at temperatures below the transforme- 
tion temperature brittle inelastic behavior is ob- 
served. This) phenomenon is exemplified in a 


remarkable fashion by sulfur vuleanizates of natural 
rubber. For this system, at room temperature, 
normal rubberlike behavior is observed when rele- 
tively small amounts of sulfur have reacted. On 
the other hand, as the amount of bound. sulfur 
increases the polymer becomes hard and inelastic. 
It is, well known that with the increase in 
the amount of combined sulfur the glass temperature, 
T,, increases from —69° to temper ee which are 
well above room temperature [1]! and it is this fact 
that causes the observed changes in siesahiea prop- 
erties 


of course 


consequence of the 
pressure is varied 
segments of 


Glass formation occurs as a 
fact that as the temperature or 
the time necessary for molecules or 
molecules to rearrange from one local liquid equi- 
librium structure to another becomes comparable 
with, or greater than, the time seale of the experi- 
mental observations [2, 3]. One important problem 
has been that of being able to qui intitative ‘ly deseribe 
the variation of 7, of a given polymer with struc- 
tural or compositional changes that are introduced. 
A great deal of success has been achieved in solving 
this general problem by considering the changes 
that occur in the free volume of the system as a 
result of the alteration of the structure or composi- 
tion [8 to 7]. Based on the theory of activated 
processes, the free volume can be considered to be a 
measure of the ease with which the necessary molec- 
ular rearrangements occur. Furthermore, the vari- 
ation of the free volume should be intimately related 


to the changes that occur in the specific volume of 


the system 


Figures in brackets it te the literature references at the end of this paper 


\ qualitative interpretation of these results is given in terms of the 


The earlier work MePherson [8S] and Seott [9] 
indicates that one consequence of the reaction of 
sulfur and natural rubber is a significant and sys- 


tematic change in the density of the svstem. These 
results give indication that the change that occurs in 
7, might be explicable by considering the variation 
of the specific volume of the system with temperature 
and composition. The data reported, however, do 
not cover a sufficiently wide range in temperature to 
be suited for this purpose. Thus there is no com- 
plate set of data available which deseribe the vari- 
ation of 7, with composition and give the specific 
volume-temperature relations for the same 
of compounds over a wide composition range 

The present work was undertaken to investigate 
this problem more thoroughly. Detailed studies of 
the volume-temperature relations of purified natural 
rubber, natural rubber-sulfur mixtures, and of vulean- 
izates of these two substances, encompassing the 
comple te ‘ly accessible composition range and cover- 
Ing #2 very wide range in temperatures, were m: nde, 
The present paper discusses the variation of 7, with 
composition in both the mixtures and vuleanizates 
and interprets it in terms of the changes that occur 
in the specific volume of the system. 


series 


2. Experimental Procedures 
2.1. Materials 


The natural rubber used in this study was a sample 
of smoked sheet, technically classified yellow. This 
sample was purified yy extracting with acetone for 
S days; it was then dried in vacuum, and 1.25 percent 
of phenyl-6-naphthvlamine was added as an anti- 
oxidant. Two master samples containing 30 percent 
and 32 percent sulfur by weight were prepared by 
mixing the rubber and sulfur on a small mill. The 
other mixtures and samples for vulcanization were 
prepared on the mill by adding an appropriate 
amount of the purified rubber to one of the master 
samples 

The vulcanization was accomplished by curing in 
a conventional press for 20 hr at 150°C. Eighteen 
vuleanizates of varving compositions were prepared, 
and the order of preparation and measurement of 


14] 








their respective physical properties was randomized 
so that any slow oxidative processes which might 
occur during storage would not bias the results. To 
further minimize oxidation and degradation all of 
the rubber samples were stored in a solid-carbon 
dioxide chest. The unreacted mixtures of the rubber 
and sulfur were prepared by pressing In a vacuum 
mold at 100° © prior to their physical examination. 
The percentage of sulfur chemically bound in the 
vuleanizates was determined by analysis after the 
removal of free sulfur by acetone extractions. The 
method of analysis was a micro Carius combustion 
followed by the gravimetric determination of the 
sulfur as barium sulfate. It was found that the 
percentage of bound sulfur was usually less than the 
sulfur contained in the mixtures prior to vuleaniza- 
tion 


2.2. Methods 


The linear thermal expansion of the mixtures and 
of the vuleanizates was determined by using a photo- 
electric recording interferometer. The design and 
operation of this instrument have been described 
by Work [10]. He concluded that 7, can be 
measured with a precision of +1.0° C with this 
instrument and that the coefficients of linear expan- 
sions are reproducible within +5 percent. 

In making these measurements three specimens 
about 5 by 5 by 3 mm were cut from each sample and 
used as spacers between the two quartz optical flats. 
The specimens were first heated in the interferometer 
to the maximum temperature to be used, then sey eral 
runs on both heating and cooling cycles were made. 
The heating was done electrically and the cooling by 
liquid nitrogen. The rate of temperature change 
was about 1° C per minute and a temperature interval 
from LSO 140° C was studied. 

The record from the photoelectric recording inter. 
ferometer gives thermocouple emf horizontally as a 
measure of temperature and a vertical step for each 
fringe or half wavelength change in thickness of the 
The linear expansion at any tempera- 
ture 7 can be found from this record by means of the 
equation: . 


to 


specimens 


ALL NX/21,+ ¢ | 

The Lomperature 7’ is determined from the thermo- 
couple emf; AL is the change in length between the 
reference temperature 7), and | i lo is the average 
thickness of the specimens at 7) (25° C in the present 
work): .V is the number of fringes between 7) and 7: 
Ais the wavelength of the light used in air at tempera- 
ture 7): and Cis a factor correcting for the change in 
wavelength mn air with temperature. The values of 
(at various temperatures for the helium vellow light 
A= 5876 A) used are given by Merritt [11] 

In interpreting the experimental results, it is more 
convenrent to consider the specific volume and its 
with temperature rather than the change in 
To accomplish this conversion the specific 


changes 

length 
Phe chet naly vas perfor! by R. A. Paulson and Miss L. J 
‘the National Burean of Standa: 
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volume v at temperature 7 was calculated from the 
linear expansion by the relation 


v=(1+AL/Ly) (2 


where 7’ 1s the specific volume at a reference tempera- 
ture 7). This relation is exact when the expansion is 
isotropic. This assumption appears be well 
founded for amorphous polymers, and the volume 
expansion calculated from interferometric data 
compares well with the direct determination of this 
quantity by dilatometric means [12]. The specific 
volumes of the samples were determined at 25° C by 
the method of hydrostatic weighing [13]. 


to 


3. Results 


The volume-temperature relations, determined in 
the described manner, for the of rubber- 
sulfur vuleanizates encompassing & composition range 
of 0 to 30° percent bound sulfur are plotted in 
figure 1. The curve pertinent to a specific vul- 
canizate is designated by a number corresponding 
to the percent bound sulfur. The glass tempera- 


series 


ture ip is taken at the pot of intersection of the 
two straight lines which represent the volume- 
temperature coefficient in the liquid and in the 
glassy state. Above 7, the volume-temperature 


data are well represented by a linear relation. How- 
ever, below 7, the linear region covers only aw very 
small temperature interval and the data are 
represented by a curve of decreasing slope with 
decreasing temperature, has been observed in 
other polymeric systems {7, 12]. Although for this 
reason a precise determination of 7, difficult 
values which are reliable on a comparative basis can 
be easily obtained. 
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Values of 7, for the various rubber-sulfur vul- 100 FT 
canizates are listed in the fifth column of table 1. 
For this series of compounds 7, continuously ‘a a 
increases as the combined sulfur content is increased, | —~ Fi | 
ranging from —69° C for the pure natural rubber to 
+-90° C for the sample containing 30 percent sulfur. | — 60 
Thus a very wide range in glass temperatures can | 
be achieved merely by reacting an appropriate / 
quantity of sulfur with natural rubber. The manner ~ 
in which 7, changes with sulfur composition is 
illustrated in figure 2. The curve in this figure is - 
distinctly S-shaped since the initial amounts of | ~ 
bound sulfur cause only a relatively slight rise in 
T,. However, when about 4 percent sulfur has | ‘; © 
been combined the rate of increase accelerates and ‘ Val 
finally levels off at the higher sulfur contents. The i 
general shape of the curve in figure 2 is similar to | ~*° 
that given by Boyer and Spencer [1] where values 
of 7, as a function of sulfur content, determined by ait 
a wide variety of methods, were collected from the ” 
then existing literature. A similar plot is also “ae 
given by Schmieder and Wolf [14] where the values -60 * a 


of 7, were determined by dynamical mechanical 
methods. In the latter work the values of 7, are 
significantly higher over the complete composition 
range than those reported here, while the compila- 
tion of Boyer and Spencer appears to be in good 
agreement with the present work except at the very 
highest sulfur contents where the literature values 
appear to be about 10° lower. 

In figure 1 the volume-temperature curves depict- 
ing each of the compounds are displaced to lower 
values as the sulfur content is increased. Thus at a 
reference temperature either in the liquid state or in 
the glassy state the specific volume ore with 
increasing sulfur content. This decrease in specific 
volume and the concomitant increase in T. | is similar 
to that which is observed in many other polymeric | 


systems [3,7]. Taking 25° C as a reference tempera- 
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FicureE 2. Plot of qlass te mperature, Ze as function of bound 


sulfur content for rubber-sulfur vulcanizates. 


ture in the liquid state, the specific volumes corre- 
sponding to this temperature are given in table 1 
and are also plotted in figure 3. For those systems 
which are actually in the glassy state at the reference 
temperature the appropriate values of the specific 
volume were obtained by a linear extrapolation of 
the liquidus. From figure 3 it can be seen that the 
specific volumes of the vuleanizates at the reference 
temperature decrease linearly with composition until 
about 19 percent sulfur at which point the rate of 
| decrease of the specific volume becomes appreciably 


TABLE 1, Thermal e2 pansion data and qlass transition te mperatures for natural rubber-sulfur vulcanizates 
Vas] a, 
Nominal Bound \ Specific I ai Slope volume 
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this change in slope would have been even larger [8] 
because of the fact that at compositions of 19 percent 
or more sulfur the glass transition is above 25° C. 
In addition to the decrease in specific volume that 
is Observed, there is also a significant continuous de- 
crease in the volume-temperature coefficient: above 
a while a definite but smaller decrease occurs in the 
linear portion below 7, The wppropriate values of 
the volume-temperature coeflicients are also elven 
in table 1 

The volume-temperature relations for some typical 
natural rubber-sulfur mixtures, as contrasted to the 
chemically combined systems, are illustrated in figure 
| For each of the mixtures the data can again be 
well represented by a straight line in the liquid re- 
gion and by a curved line in the glassy state. The 
glass temperatures, as determined from the data of 
figure 4, are remarkably constant and can be assigned 
a value of —68° C +1° for all the mixtures. The 
dependence of the specific volume of the mixtures at 
Zo” 6 © 


on the sulfur content is also shown in figure 3; 
The 
values of the specific volume observed and those 
calculated by the additive mixing law agree very well 
and indicate that no net change in volume occurs on 
mixing. A similar observation was made some vears 
avo by MePherson [S|]. The pertinent physical con- 
stants for these mixtures are tabulated in table 2 


the data are well represented by 2 straight line 


4. Discussion 


The very large changes in the glass temperature of 
natural rubber that are induced by with 
sulfur must be caused by the alteration of the liquid 
structure 


reaction 


State structure These changes ll occur 


as al conseqtuence of the chemical processes that have 
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taken place. The natural rubber-sulfur reaction, 
though well known and intensively studied, is quite 
complex. There are wide areas of disagreement 
among various investigators as to the mechanism of 
the reaction and the resulting structures that are 
formed [{15, 16]. It is clear, however, that a portion 
of the combined sulfur is utilized to form intermo- 
lecular crosslinks. The complexity of the chemical 
precludes at present the possibility — of 
ascribing the changes in 7, to specific chemical 
alteration, and somewhat more general considerations 
must be given to explain the drastic changes in 
the physical and mechanical properties that occur. 

The variation of the glass temperature in other 
polymeric systems with the svstematic alteration of 
structure or the composition has been satisfactorily 
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TABLE 2 Thermal expansion data and ala fransitior 
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explained on the basis of free volume considerations 
'3-7|. The problems treated, using this concept, | 
have been those where the free volume itself need 
not be specified but the changes in the free volume, 
as manifested by the changes in the specific volume, 
are determined, and the change in 7, relative to that 
for a reference polymer can be calculated. This 
procedure is limited to cases where the changes that 
are occurring do not drastically alter the chemical 
or molecular structure of the polymer and has been 
successful in treating the molecular weight depend- 
ence of 7, |3], the effect of crosslinking [4], dilution 
with low-molecular-weight materials [5], copolymeri- 
zation |5, 6], and the effect of increasing the number 
of groups in a polymer side chain [7]. It is of interest 
to apply the aforementioned ideas to vulceanizates 
of rubber-sulfur. 

In utilizing this method the specific volume in the 
liquid state is assumed to consist of two parts; that 
volume which would be occupied by a closely packe dl 
liquid structure and the remaining volume which is 
termed the free volume. The net change in the spe- 
cific volume of a given system, when subject to an 
external change, may reflect changes in the occupied 
volume, in the free volume, or in both. One conse- 
quence of the reaction of natural rubber with sulfur 
is that at a reference temperature in the liquid state 
the specific volume is reduced. Following proce- 
dures used in deseribing 7, for the poly(n-alkyl 
methacrvlate) series of pols mers, a certain fraction / 
of this volume change is attributed to a decrease in) 
the free volume of it svstem, 
change in free volume that occurs at temperature 7 
can be expresse “das 


T) —fle°(T) —3"(T)] (3) 


where (7) and ¢°(7) are the specific free volumes 
of the vulcanizate and pure rubber, respectively, and 


B°(T) and 7°(7) are the corresponding values of the 
specific ot tink A at a temperature 7 in the liquid 
state. It is further assumed that the free volume 


varies linearly with temperature so that 
¢o( 1) =¢,+ Aa( T— T,) (4 


where ¢, is the free volume at the glass temperature 
7T,, and Aq@ represents the rate of change of the free 
volume with temperature. By writing equation (4) 
for the vulcanizate and for pure rubber at a tempera- 
ture 7’, and substituting these results in equation (3), 
the following expression is obtained relating the 
glass temperature of the vuleanizate 7 to that of 
the pure polymer and to the volumes of the system: 
Te = Tp—(Aa®/Aa‘)(Tp— T?) 


t 


(f/da')[B%T) 
ne Tp)| ‘ (@,° od, ) Aa (5) 
where 7, is a reference temperature taken to be 


298° K and the superscript ° refers to pure rubber 
while the superscript * refers to the vulcanizate. 


so that the relative 


| 


n order for equation (5) to be converted into a 
mm... that can be useful in interpreting the experi- 
mental results, it is necessary that the appropriate 
criterion for formation be specified. The 
simplest condition to take is that glass formation 
represents a state of isofree volume. This criterion 
was originally applied by Fox and Flory [3] in 
explaining the effect of molecular weight on the glass 
temperature of polystyrene and has been applied 
in explaining the effect of composition in copolymers 
and polymer-diluent mixtures and the effect of 
structural changes on the glass temperature of the 
poly(alkvl methacrvlates). When the aforementioned 
criterion is applied, equation (5) simplifies to 


glass 


Ty =298— 94 Aa’ / Aa’ + f/ Aa‘[v"(298) — v*(298)] (3) 


It has been suggested [17, 18] that the condition 
o,/?, Is a constant at 7, is a more. satisfactory 
specification of glass formation. When this condi- 
tion is introduced into equation (5) an additional 
term results in equation (6). However, it has been 
found that this slight correction does not sensibly 
alter the results |7]. 

Utilizing equation (6), values of f for each of the 
vuleanizates can be computed by using smoothed 
data derived from table 1 and assuming that the 
volume-temperature coefficient for the free volume 
can be identified with the difference in the cor- 
responding coefficients in the liqnuid and 
state. The quantities used in the calculation and 
the resulting values of f that are computed from 
equation (O) are viven in table 3. It can be seen 
from these results that a reasonable value for the 
magnitude of the quantity f is obtained; it- being 
less than unity and comparable to that obtained 
for other systems. As contrasted to the results for 
other svstems, however, these deduced values of f 


olassy 


TABLI 3 (Calculation oft f. the fraction. of the speci fu Volidiie 


chanue upon vulcanization which occurs in the free olume 

Bou = - 

ulfur x JUS JUS 
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Che values in the second and third columns are taken from smoothed curves, 
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remain only approximately constant with increasing 
amounts of bound sulfur. They gradually increase 
from a value of 0.23 at 2 percent combined sulfur 
to a maximum value of 0.32 at 17 percent sulfur 
and then slowly decrease once again. Though the 
deviation of each value of f from the mean is not 
very large, the lack of constancy in this quantity is 
such that using the mean value of f in equation (6) 
does not adequately represent the data of figure 2. 
From free volume considerations it would be ex- 
pected that the glass transition temperature for the 
rubber-sulfur vulcanizates should increase the 
amount of combined sulfur increases since a concomi- 
tant decrease in specific volume occurs. Though this 
effect is observed, the simple relation that has been 
derived in equation (6), though known to be appli- 
cable to other systems [7], is only qualitatively cor- 
rect for the system discussed in the present work. 
The reason for the lack of quantitative agreement in 
this instance, when compared to other systems, can 
be due to complications caused by structural altera- 
tions induced by the chemical processes. Even at 
the lower levels of bound sulfur only a portion of the 
sulfur involved in intermolecular crosslinks, the 
fate of the remainder of the sulfur being still a matter 
of argument [16]. The structure of the compounds 
containing the higher amounts of bound sulfur is still 
With the distinet possibility existing of 
major structural changes on reaction with sulfur it is 
not surprising that the relatively simple considera- 
tion set forth above vields only qualitative agreement 
with experiment. It must also be remembered that 
$2 percent bound sulfur corresponds to one sulfur 
atom bound per isoprene unit. On this basis a rela- 
tively large amount of sulfur is incorporated into the 
compounds and could thus be treated as a second 
component, similar to a copolymeric ingredient or 
low molecular weight diluent. Although, in a formal 
fashion, the problem can be treated from this point 
of view by methods previously (5, 6], 


as 


Is 


obscure. 


deseribed 


several arbitrarily selected parameters are required 
which can not be evaluated independently. Thus, 
the physical significance of the results is open to 
question. 

On the other hand, in the case of the natural 
rubber-sulfur mixtures the concepts outlined above 
appear to apply very well. In this instance there is 
no net change in the specific volume of the system on 
mixing over the complete composition range. Thus 
the free volume available to the polymer is neither 
increased nor decreased and consequently 7, remains 
invariant with composition. 
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ral Elastic Deformations in Strips With Holes Loaded 
is Through Pins 
on 
ius 
ler Michael Chi and L. K. Irwin 
Ins 
Bickley’s solution for the generalized plane stress in a plate loaded through a pin is ex- 
tended to give the deformation at any point in a plate loaded through one or two pins which 
are on the same load line. Analytical expressions are derived that are known to apply to 
a finite strip with a width >14a and distance between holes >a, where a is the radius of 
id the hole. Numerical examples are developed and presented in the form of graphs for typi- 
2 cal cases. 
6). 
”). 1. Nomenclature 2. Introduction 
J a=radius of the hole The loads applied to a plate can be transmitted to 
K,,=constant coefficients, e.g., Ky=', Ay ‘2, | other plates through joints which are made with 
1, k,=0, K3=20, As=0, ete. fasteners in holes. The loads produce shearing 
E=Young’s modulus of elasticity | stresses in the fasteners which, in turn, produce com- 
0; G=modulus of rigidity, ae | pressive or “bearing” forces on the boundaries of the 
~? 2(1+p) | holes. The elastic stresses in the plate produced by 
n=positive integers from 0 to © in an infinite | the boundary force on a hole, if solved, can be used 
J. series whose general term is /(7) | directly for the determination of the elastic deforma- 
n,—a positive integer arbitrarily chosen to desig- | tions. These deformations are of primary importance 
“h nate the first term of a partial series | in evaluating the suitability of different Combinations 
| p—maximum intensity of the pressure from the | of materials and fasteners for aeronautical structures, 
~ pin on the hole | This method cannot be applied to a strip of arbi- 
“e P=total force exerted by pin on the boundary | trary dimensions loaded by pins because the general 
of hole | solution for the elastic stresses in this case has not 
: r.@ polar coordinates | been found. However, the problem of the elastic 
, s=partial sum of a series stresses in an infinite plate loaded by one rigid pin 
u—displacement in radial direction was solved by Bickley [1]' and the problem of a 
r=displacement in tangential direction plate with a circular hole at its middle subjected to 
w,—displacement in longitudinal direction | uniform unidirectional tension was solved by Kirsch 
wy, —displacement in transverse direction (2). The stresses given by the solutions contained in 
Yro—shearing strain these papers are combined by superposition to ob- 
€-—strain in radial direction tain the stresses for the problem considered here. 
€9—strain in tangential direction With the aid of these stresses, the deformations of 
u— Poisson’s ratio finite strips loaded through one or two holes located 


u’ =modified Poisson’s ratio in a generalized plane | on the load line are caleulated. 

stress problem, u/(1+ 4) 
o,—normal stress in radial direction 3. Theory 
og—Nnormal stress in tangential direction 
T@— shearing stress 





Bickley’s solution for the general plane stress func- 
| tion @ for an infinite plate loaded by one solid pin is: 


. F , r , l r r 
o a°| As log, + B,6+-A, { ésin é@ = ( log, -) eos o} 
a a 2a - 


r l Fy r : , a cK 
+B ¥! #@ cos 6 — (log, )sin ob 5 cos 6 Dy sin @ 
a o a a r , 


+S fart promt (4, cos mé+ B,, sin mé) +a" r-” (C,, cos mé+-D,, sin mé) ] (1) 
’ . Figures in brackets indicate the literature references at the end of this paper 
where «1,,, B,,, ©, and D, are coefficients. | se , 
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The above expression satisfies the biharmonie | half the boundary of the hole [1; also 4, p. 123]. This 
equation based on compatibility considerations [3, | has been confirmed by experiment [1]. Since cosine 
p. 54], functions are even functions, we can represent the 

distribution of the boundary forces on the edge of the 
Oo 10,1 0&\/0%@,10¢, 1 O%@ hole due to the pin by the Fourier cosine series: 
or or’ Oe )\ or? ror) 2 Oe? ) 
, :, f(0) — "4+ Sra, cos mé. 3 
and has single-valued deformations and stresses » il, a 

The stresses in polar coordinates in terms of the 

stress function @ which satisfy the equations of | The coefficients are determined by a known distribu- 





equilibrium are [3, p. 52]: tion of boundary forces, 
ldo, 1 0°) T T 
0 t - : 1(@ p cos ff. <~.7— 
rar fF Ce 2 9 
O*¢ 0 lsewhere 
1 ? ‘ 9 elsewhere 
O/ 
Since the stresses at the edge of the hole must be 
O | oF equal to the boundary forces, due to the pin, the co- 
— or\r Ob/ J efficients in eq (1) may be determined by substituting 


eq (1) in eq (2), performing the indicated partial dif- 
ferentiations, and then equating the resulting stresses 
at the edge of the hole to the forces, eq (3). If the 
distribution, /(@), is assumed acting, the stresses, eq 


The distribution of the forces due to the pin load is 
customarily assumed to Vary as a cosine function ovet 


2 . become 











” ‘5 
| (s ) »> ua | rv a 
( } + Cos cS 
ah } Ss / s & ) ” 
n | . mm | 
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Since eq 1) was derived using plane strain theory 
the above equations were made applicable for a gen Ou | » 
eralized plane stress solution by replacing Poisson’s 0 € ke 0 uo 7 
ratio, wu, with w’ which is equal to w (1m) [5 , 
In pole coordinates, the stress strain deformation O/ / 
relationships ure On re r€e Ke On ua > ) 
lou oOo | 
r Ov O/ / ; (y J 
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Substituting the first and second expressions of (4) 
into the first part of eq (5) and integrating, 


u | ca al mI \d, ue’ oe) dr 


similarly, 
a | ; {7 fi : lr 6 
/ | / 0a Mo,ja | \0, M G@) al (0) 
By substituting (6) into the last of eq (5), it can be 
shown that the constants of integration are rigid body 
displacements |[3, p. 189] and therefore may be set 
equal to zero for the purposes of this investigation. 
Thus far, only the infinite plate with the load ap- 
plied through a pin at the middle and supported at 
infinity [6] as shown in figure 1(a) has been con- 
sidered. The complete solution of the problem, fig- 
ure l(c), is obtained by superposing the solutions of 
two subproblems, figures 1(a) and 1(b). The prob- 
lem of a plate with a circular hole at its middle and 
subject to a uniform unidirectional tension, figure 
1(b), is the well known Kirsch’s problem which can 
be found in standard texts [3, p. 195]. 


pee heeaabaaa aaehaebhabanaa eaebbhbadbbadr 


——-— - a. a | 


——.— 4 
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FIGURE 1. Superposition of subproblems (a) and (b) to solve 
proble ne (CG). 


Sketches represent infinite plat 


Combining these two solutions, the equations for 
radial and tangential displacements are, respectively, 


- ; ' ‘ ; 9 ‘\9 ‘ 5 
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w , 


is taken 


as 0.25 for the purpose of this investi- 











gation. 

The displacements in Cartesian coordinates can 
be evaluated from the radial and tangential displace- 
ments in polar cordinates by the following relation- 


ship [4, p. 18]: 


w, “wu cos 6—v# sin @ 
\) 
w usin 6+” cos @. 


Therefore, the displacements at any point in the 
infinite plate loaded through a pin can be found by 
substituting eq (7) into eq (8). 

The displacements thus obtained are meaningless 
by themselves because they include an arbitrary 
amount of rigid body displacements. The displace- 
ments of interest are the relative displacements be- 
tween the points in the plate and a chosen reference 
point. The point, r=a, @=0 is chosen as the 
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reference point. The displacement of any point Is 
taken to be the difference between the displacement 
of that point as computed by eq (8) and the dis- 
placement of the reference point. 

The rate of convergence of the series in the equa- 
tions for the displacements, eq (7), varies with the 
distance of the point from the hole. The first five 
terms were sufficient to give an accuracy within 
one-half percent for the displacements at all points. 
The relative displacements are the differences be- 
tween large numbers of the same order of magnitude. 
To insure sufficient accuracy in the final results, the 
displacements in eq (7) were therefore carried to 
nine significant figures. This involved, for example, 
the inclusion of 6 terms of the series for the points 
a distance 3 times the radius of the hole, 3a, from 
the edge of the hole and 16 terms for points a dis- 
tance 0.5a from the edge of the hole. For comparable 
accuracy at the edge of the hole, about 100 terms of 
the series would be required. To reduce the com- 
putational task, the following Euler-Maclaurin For- 
mula for computing the partial sum of the series was 
employed. 


" ] 
i ®t] (d” : ’ 
s=S* fin) 4 f(n)dn +>) 4,2” a J n) | 
ale . ; m= adi ' 


2 4n-+-l—p 


In this case, f(n) re 


wr (40° 


4. Finite Strip With One Hole 


In a plate of finite width, a hole close to an edge 
of the plate may have considerable effect on the 
stress distribution in the plate. The case of a plate 
of finite width (strip) with a circular hole on the axis 
of a symmetry was solved by a method of successive 
approximations [7]. The method was later extended 
to the case in which the hole is loaded through a 
pin [8S]. However, the strain expressions obtained 
by this method were not suitable for integrating to 
obtain the displacements. 

For strips whose width is large compared with the 
diameter of the hole, ratio of diameter of hole to 
width of strip less than 0.15, the effects of the edges 
of the strip on the stresses around the hole can be 
neglected [9]. On this basis, the displacements of an 
infinite plate due to a load applied through a pin are 
assumed to be applicable to a finite strip without 
serious error Figure 2 shows the relative displaces 
ments along the longit udinal centerline and the edge- 
for a strip with one pin. The diameter to width 
ratio is equal to 0.125 


5. Finite Strip With Two Holes 


The stresses around a hole are influenced by the 
presence of adjacent holes. The method of successive 
approximations is again applicable to this case [10] 
but becomes much more complicated for the case of 
two or more holes with boundaries loaded through 
pins 
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Figure 2. Relative displacements along the longitudinal center- 


line and the edges for a strip with One pin. 


As previously mentioned, the effect of a hole in a 
plate is a local phenomenon, and the effects of ad- 
jacent holes on each other can be neglected [9] if 
they are sufficiently far apart; that is, if the center to 
center distance of the holes is four or more times the 
diameter of the larger hole. Thus for a spacing 
exceeding four hole diameters, the displacement of a 
point in a plate containing two holes loaded through 
pins can be taken equal to the sum of the displace- 
ments due to each hole taken separately. 

The relative displacements along the longitudinal 
centerline and along the edges of a strip with two 
equidiameter holes located on the longitudinal center- 
line and loaded through pins are shown in figure 3 
The holes are four diameters apart and the ratio of 
the diameter of the hole to the width of the strip 
equals 0.125. 


6. Concluding Remarks 


The analytical expressions and the computational 
methods presented here are applicable for determin- 
ing the elastic deformations of finite strips loaded by 
pins when the effects of edges and interference be- 
tween adjacent holes are negligible. Such strips 
have widths > 14a and distances between holes > Si. 
The examples for which deformations are shown in 
figures 2 and 3 correspond directly to single-shear test 
specimens widely used to evaluate the performance 
of joints. Further studies should be made to extend 
the range of these results. 


This investigation was part of a project conducted 
under the sponsorship of the Bureau of Aeronautics. 
The authors wish to thank this agency for its support 
of this work. 
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Absolute Light-Scattering Photometer: 
I. Design and Operation 
Donald McIntyre and G. C. Doderer ' 


\ new light-scattering photometer has been designed and built for determining the 


absolute seattering from polymer solutions. 


The instrument is also capable of performing 


as a research instrument for making measurements at very low and very high angles, and 


at very low and very high intensities of scattered light. 


The instrument seans the angular 


scattering either manually or automatically while measuring continuously the ratio of the 


seattered light to the incident light. 


1. Introduction 


The general need for accurate light-scattering 
measurements for the determination of molecular 
weights of high polymeric compounds is well known.* 
Some of the problems connected with these measure- 
ments and most of the work up until 2 vears ago 
have been summarized in two publications.’ * The 
present work was done with the object of increasing 
the accuracy of determining molecular weights of 
macromolecules rather than that of determining 
optical constants of materials. 

The instrument deseribed here was designed to 
attain an accuracy of 3 to 5 percent in the determina- 
tion of molecular weights. An instrument that must 
be versatile enough to make molecular weight meas- 
urements in any research problem and also make 
absolute scattering determinations imposes many 
compromises in design. For example, a closed-in 
system is much more convenient for general labora- 
tory work, whereas an open system in a darkroom 
would be more suitable for absolute measurements. 
A compromise is to make a rather large working 
compartment in the instrument. Also, an absolute 
measurement would only have to be carried out at 
a few fixed angles so that considerably more rigidity 
and inflexibility could be built into an instrument 
to attain more accurate measurements. In con- 
trast, a research instrument must be capable of 
scanning a large angular range with a measuring 
technique that is not unduly time consuming and 
fatiguing. 

With these thoughts in mind an instrument was 
constructed with emphasis on providing (1) a variety 
of precise experimental setups involving different 
combinations of sources, lightbeams, filters, sample 
cells, and receiver arrangements; (2) a well-colli- 
mated beam of monochromatic light having uniform 
intensity and producing negligible stray light inside 
the instrument; (3) a sensitive and stable electronic 
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detecting and recording system; and (4) a rugged 
mechanical design featuring easily adjusted com- 
ponents. 


2. Apparatus 
2.1. General Description 


The over-all plan of the instrument is shown in 
figure 1. A transformer connected to a voltage 
regulator operates a mercury vapor lamp to produce 
a small intense source of light. The image of the 
source is then projected as a monochromatic beam 
of variable collimation that has well-defined dimen- 
sions at the center of the cell (C). From there the 
beam passes on to a light trap where it is absorbed. 

A very small percentage of the projected beam is 
reflected by a cemented glass prism on to a monitor- 
ing photomultiplier tube (Mon PM). The measur- 
ing photomultiplier tube (Meas PM) is mounted 
on a T-beam connected to a goniometer so that it 
can be rotated to view the light scattered from the 
sample, which is mounted on a stationary pedestal 
at the center of rotation. Both phototubes are 
similar and operate from the same power supply 
controlled by a voltage regulator. The phototube 
outputs are amplified and fed to a ratio recorder 
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which is essentially an automatic precision potenti- | Corning 4600 special heat absorbing filter F;. 


ometer operating a chart recorder. The recorder 
compares the voltage from the measuring circuit with 
the reference voltage from the monitoring circuit. 
Thus the amount of scattered light is measured 
relative to the incident beam and ts recorded as a 
value which is independent of the absolute beam 
intensity. 

Figure 2, an exterior view of the instrument, shows 
the source enclosure, optical projector, and goni- 
ometer housing mounted on a table and connected 
to the electronic control panel. 





FIGURE 2. 


Overall view of phe tometer 


2.2. Optical System 


The essential optical components are represented 
in figure 3. Two source lamps have been tested: an 
S5-w mercury vapor lamp, and a 250-w compact- 
source medium-pressure mercury vapor discharge 
lamp. The latter forms a discharge approximately 
$3.75 mm high by 1.5 mm wide between two chisel- 
edge electrodes. The chisel edges are alined in the 
direction of the optic axis so that the are wanderings 
tend to occur along, rather than across the axis 

The source is located at the focal point of a plano- 
convex condensing lens L, of 84.5 mm focal length 
which receives light through a wide stop S,; and a 
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Essen- 
tially parallel light from lens L, passes through 
removable filters F, and F;. F, represents three 
monochromatic filter combinations: (1) Corning 
glass filters 3389 and 5113 passing mercury blue light 
of wavelength 436 mu; (2) Corning glass filters 3484, 
5120, and 4303 passing mercury green light of wave- 
length 546 my; and (3) a multilayer interference filter 
also passing the mercury green line. F; represents 
1 mm thick optical neutral glass filters, as do Fy, 
and Fs. 

A second condensing lens L., similar to L,, focuses 
an image of the source on aperture A, located at the 
focal point of lines Ls. The aperture A, is fitted into 
a milled slot and is removable so that variable degrees 
of beam collimation are obtainable. The largest 
aperture is 3.4 mm high by 1 mm wide so that it ts 
smaller than the source image. Field stops which 
determine the size and shape of the projected beam 
ia be inserted at position S. and adjusted axially so 
that the achromatic collimating lens L, of 98.6 mm 
focal length and 40-mm diameter focuses a sharp 
image of the field stop in the center of the sample cell, 
C. The beam can be shut off at the entrance to the 
main instrument housing by shutter Sh,. The rela- 
tive positions of the above components may be seen 
in figure 4. The remaining optical components in- 
side the main housing are shown also in figure 5. 

A l-in. glass square consisting of two 45°-90 
prisms cemented together with Canada balsam is 
indicated as G in figure 3. The cemented interface 
reflects a very small percentage of the incident light 
into the monitoring system. By choosing glass and 
cement of nearly the same refractive index, the same 
beam as that viewed by the measuring phototube 
may be monitored with minimum intensity loss and 
negligible polarization. Using barium crown glass 
of mg 1.5415, which was obtained from samples made 
at the NBS glass plant, and Canada balsam, which 
was reported to have an ng— 1.5446, the measured 
reflectance from the surface was 10~° of the incident 
green light. Thus the match of the refractive indices 
at this wavelength is closer than the above figures 
would indicate 
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FIGURE 4. Detailed view of the projector. 


The main beam is attenuated by movable filters 
F, having transmittances of approximately 1/2, 1/3, 
1/10, 1/40, and 1/4,500. P, and P, are Polaroid 
KN36 glass-laminated filters mounted with axes 
horizontal and vertical, respectively, and having 
maximum deviation of 4 min and maximum resolving 
power of 1 min. Shutter Sh, allows the measuring 
beam to be shut off without disrupting the moni- 
toring system. S; is a baffle to eliminate stray 
light. In addition, removable baffles can be inserted 
into this nosepiece arrangement as beam dimensions 
are changed. Another tube can also be inserted 
on the cell table which is capable of accommodating 
two additional baffles. Thus there are four possible 
interchangeable baffles in addition to two shutter 
diaphragms, which are circular for eliminating stray 
light. A 4-mm-thick neutral glass filter Fy of approx- 
imately 1/30,000 transmittance can be mounted on 
the end of the measuring phototube arm and subse- 
quently brought into the lightbeam when the photo- 
tube is rotated into the beam. This is mounted as 
a milled slot so that other neutral filters canbe 
inserted as they are needed. This arrangement is 
used only for automatic scanning of the entire 
angular scattering, and for special purposes of 
calibration 

The receiver optics consist of an achromatic lens 
L. of 40.5-mm focal length and 18.5-mm diameter, 
a small aperture A, mounted at the focal plane, and 
a stop S; which determines the area of lightbeam 
that is to be admitted to the phototube PNI,. 
There is a series of aperture stops, Ag, that can be 
easily and quickly inserted to define the angular 
acceptance of the receiver, just as there is a series 
of field Stops, Os, that can be inserted into the 
receiver in order to change the field of view consonant 
with the physical limitations of the phototube. 
The position of the whole receiver is adjustable so 
that the image of the field stop can be placed directly 
at the center of rotation. A Polaroid HN22 filter 





Detailed view of receivers and cell table 
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P; is mounted in a slide and can be adjusted to 
rotate 90° between the axes of the polarizers. A 
second slide I. may be used to hold a fluorescence 
filter. In addition, provision is made to insert a 
diffusor behind the field stop when necessary, 

The monitoring optics consist of lens Ly, aperture 
Ao, field stop Sy, and phototube PM,, similar to those 
in the measuring receiver just described. Attenuat- 
ing filters Fy having transmittances of approximately 
4% and !%5 are also placed before this receiver. 

There are two interchangeable phototube housings 
on the measuring photo tube. Each of these housings 
can accommodate the receiver mount described, and 
placed into position without significant 
changes by means of milled slots. One housing con- 
tains a 1P21 photomultiplier which is used routinely, 
while the other housing contains an end-on photo- 
multiplier, 5819, for special work. In addition, the 
latter housing can take a cassette which uses 1 by 3 in. 
photographic plates in place of the photomultiplier. 


can be 


2.3. Electronic System 


The sensitive and stable detecting-recording syvs- 
tem is based upon the use of matched precision com- 
ponents in the monitoring and measuring circuits 
which operate the ratio recorder. Optionally, the 
recorder may be operated as a typical potentiometer 
recorder to indicate the output of either phototube 
separately . The receivers are closed desiccated metal 
housings containing RCA 1P21 multiplier phototubes, 
A power supply having l-ma d-c output, stabilized 
to 0.5 percent and fed by a 0.250-Kva voltage stabi- 
lizer, is used for the high-voltage supply. When 
operated at S00 v d-c the phototube dark currents 
are about 310°" amps. The monitor light beam 
may be adjusted to its maximum unfiltered intensity, 
which produces a total output of over 1 107° amps. 
Thus dark-current fluctuations have a negligible 
effect on the final monitor reference voltage. In 
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practice, a monitor output of about 1 107° amp has 
proven very satisfactory. The phototube outputs 
are amplified by two Leeds and Northrup 9836A 
stabilized d-c micromicroampere indicating amplifiers 
providing a full scale output of 10 mv over a range of 
1107" to 210°° amp with an accuracy of +0.5 
percent of the range or 0.5107" amp at the most 
sensitive region. The measuring circuit: amplifier 
supplies voltage to the potentiometer of a Speedomax 
ratio recorder where it is compared with the monitor 
reference voltage to an accuracy of 0.3 percent of the 
range. 

Thus fluctuations in source intensity affect both 
input voltages proportionately, and the final recorded 
value is independent of the absolute beam intensity. 
In operation, the measured values are corrected by 
subtracting a dark current reading with the shutter 
Sh. closed. This is done either arithmetically or by 
the use of a biasing voltage on the measuring photo- 
tube signal using a mercury cell. 

The 85-w lamp operates from a transformer on 
110 v a-c; whereas, the 250-w lamp operates from 
a ballast on a 220 v a-c line. 


2.4. Mechanical Elements 


The instrument was constructed as three inde- 
pendently adjustable assemblies; namely, the source, 
the optical projector, and a light-tight measuring 
compartment mounted on a heavy welded 
framework table top. All of these pieces 
optically alined using the center of rotation of the 
measuring phototube and goniometer as a fiducial 
mark. Visual observations of the beam by means of 
targets placed in indexed positions along the rotating 
beam and center of rotation were used for the final 
alinement. Although the overall construction ts 
sturdy to minimize possible vibration or movement, 
a wide range of adjustability is maintained. Thus 
the receiver apertures and field thread into 
place against positioning rings, and the mounts are 
positioned with the aid of pins. 

The measuring phototube housing is mounted in 
machined ways on an aluminum beam of large T- 
With a narrow slit in front of the 
receiver lens, the entire housing Trha\ be moved across 
the lightbeam bv a screw to determine the uniformity 
of the light flux. The aluminum beam ts connected 
to a vronlometer which is eraduated in degrees from 

LO” to 150° with respect to the imeident light 
The assembly may be rotated manually 
or automatically in either direction by 
geared motor drive. A scanning rate of 40 
used to S\ nehronize with the chart recorder speed of 
2 ingmin. Two electrical limit switches prevent 
mechanical damage due to overscanning and a third 
microswitch is operated by pins spaced 10° apart on 
the goniometer to activate a magnetic pen marker 
which produces a fiducial mark on the chart. 

Sample cells and thermostats are positioned accu- 
rately on astationary pedest: alat the centerofrotation 
by a three-point support system. The pedestal top 
three ‘-in.-diam ball bearings on *;.-in.-diam 


steel 


were 


stops 


cross section 


as desired, 
a reversible 
min was 


has 


holes spaced 120° apart on a 3-in.-diam circle centered 
on the axis of rotation. Each cell assembly is cement- 
ed to a similarly drilled plate which automatically 
indexes on the pedestal with no undesirable mechan- 
ical tight fits or fastenings. 

The light trap has two sections containing dark 
neutral filter glasses mounted at 45° to the incident 
beam. These filters at this angle of incidence absorb 
most of the light, reflecting only a few percent from 
their surfaces. In the large first section which 
absorbs low intensity scattered light surrounding the 
main beam, the filter surfaces reflect any light onto 
the black-painted trap walls where further reflection 
is negligible. The center of the first section is open 
and passes the main beam into the second section 

In this section the light undergoes two absorptions 
and 45° reflections by filter glass until it is finally 
directed to a black-painted portion of the trap not 
visible to the interior of the instrument compart- 
ment. The entire light trap is removable so that the 
interior may be viewed by a camera or telescope. 

The photometer is also capable of operating at 
high temperatures since the cell table fixtures are 
located only by means of ball bearings which par- 
tially insulate the shaft below from the cell. Also 
the individual photometer housings are water cooled 
on the front plate to prevent any increase in the dark 
current by radiation from the sample. 


3. Performance 


: 3.1. Electronic 

The 250-w lamp provides a light intensity approxi- 
mately 11 times that of the 85-w lamp, but the are 
is relatively unstable and would produce unacceptable 
errors in a single channel instrument, as is shown in 
the upper portion of the recorder trace in figure 6 
The lower portion shows a trace under identical con- 
ditions, but using the monitored ratio-recording sys- 
tem. Here the trace variations are reduced to the 
small noise inherent in a high level d-e amplification 
system. From such a trace a constant average value 
is easily read. 

In order to determine the response characteristics 
of the measuring system over a large range of intensi- 
ties the transmittance of a neutral filter of approxi- 
mately ‘ transmittance was determined while re- 
ducing the lightbeam intensity with a neutral filter 
In this work a large viewing field was used on the 
photomultiplier to eliminate the small inevitable local 
variations in a neutral filter. In changing the inten- 
sity by a factor of 1,000, the measured transmittance 
of the filter showed a maximum deviation of less than 
| percent with an average deviation of less than 0.5 
percent. This is a severe test because in the very low 
intensity measurements the signal is much less than 
that from pure liquids and the ratio of dark current 
to total current becomes much larger. 

Over several hours the system is very stable; how- 
ever, the ratio of scattered light from an object will 
eventually change due to a small change of the tem- 
perature of the cemented prism. The small change 
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trace of scattered light from a diffuser with (1) a single measuring photomultiplier, and —2) a ratio of two 


photomultiplier signals 


in temperature in the prism can affect the percentage 
of reflected light so that the ratio can change. This 
is not critical since all measurements are made in 
the course of a few hours, and those kinetic measure- 
ments which last longer can easily be referred to the 
direct beam whether the scanning is manual or auto- 
matic. 
3.2. Photometric 


Stray light. In the early experiments observa- 
tions were made to determine the existence of stray 
light which might be caused by some unforeseen 
peculiarity of the apparatus. Thus experiments 
were conducted with a nosepiece containing two 
stops which allowed the phototube to view a 6- by 
10-mm field while accepting about 2.5°. The pro- 
jected beam was about 10 mm wide and 12 mm high 
with a divergence of 0.5° in the horizontal direction 
and 1° in the vertical direction. Measurements were 
initially made shortly after changing stops with the 
box open. The photomultiplier current is measured 
as a function of viewing angle and of time after clos- 
ing the top of the instrument. The data are listed 
in table | Below 12° the phototube receiver by 
visual observation appears to enter the edge of the 
large beam. The current in excess of the dark cur- 
rent can be attributed to stray light. From these 
results it appears that unless measures are taken to 
allow all of the dust to settle in the box the results 
in the forward direction will always contain a large 
amount of scattered light. 


Pasie | Photomultiplie irrents at different angles and 


fimes 


The effects of adding diaphragms to eliminate 
stray light were also determined. There are four re- 
movable stray light baffles for a given beam size. 
They are located about 6 in. away from the pro- 
jection lens and follow two circular iris diaphragms. 
They can be inserted about every 25 mm thereafter 
in the projected beam up to 100 mm from the center 
of rotation. They are in all cases slightly larger 
than the beam dimensions calculated for that loca- 
tion in the box. With the large beam (10 mm wide) 
no difference could be detected from 13° just out- 
side of the large beam to 150° away from the direc- 
tion of the beam when the stray light baffles were 
inserted or removed, Later a large cylindrical cell 
with plane entrance and exit windows and _ filled 
with a liquid whose refractive index was not too dif- 
ferent from that of the glass was placed on the cell 
table and again there were no differences when the 
baffles were in or out of position or over the above 
angles. However, when the narrow beam (3 mm) 
Was examined without a cell in place the Stray light 
decreased by about 25 percent below 10° when ap- 
propriate baffles were introduced into the system, 

In order to investigate further the effect of stray 
light the receiver was removed from the phototube 
housing so that the photomultiplier was receiving 
light from any angle that its photosensitive surface 
would allow. Again the addition of the baffles 
with the large beam while viewing from 45° to 150 
showed no appreciable change in the response of 
the phototube. However, visual inspection after 
long dark adaptation indicated that a baille in one 
of the middle positions did decrease some of the 
diffuse illumination in the photometer. From the 
results of all these experiments it is felt that this 
photometer can determine absolute amounts of 
scattered light from the usual organic liquids without 
any appreciable error from stray light 

Fluorescence, Polarization, and Monochromatic 
Filters. The directions of polarization were seribed 
on the Polaroid filters that were subsequently used 
for the polarizers. They were then mounted in 
the instrument by alining the inscribed marks with 
respect to the perpendicular and horizontal lines at 
the cell table. The analyzer was then set for 
extinction of the polarized beam. The percent 


157 








of transmission at extinction for the analyzer was 
0.012 and 0.008 for blue light (436 my) in the 
horizontal and vertical direction, respectively, and 
0.008 for 546-my light in both the horizontal and the 
vertical directions. The phototube was found to be 
equally sensitive to both horizontally and vertically 
prolonged light of 436 my wavelength, but was 0.5 
percent more sensitive to horizontally prolonged light 
of 546 my wavelength. 

Four filters for use with fluorescent solutions can 
be mounted in the barrel of the receiver. The No. 
3384 Corning filter transmitted 9107* percent 
of the 436-my light, and 8S percent of the 546-my 
light. The No. 2424 Corning filter transmitted 
1.61077 percent at 546 my and 7.8107° at 436 
Mu These are simply cutoff filters to detect the 
amount of fluorescence. A of monochromatic 
filters similar to those described below could 
be used to allow light-scattering measurements to 
be made on fluorescent solutions. 

The monochromatic filter combinations used for 
the incident beam were examined on a Beckman DK 
spectrophotometer and found to give less than 0.5 
percent of the maximum intensity when measured 
at 200 A from 4358 and 5461 A. The transmission 
of the interference filter was about 2.5 times that 
of the glass filter. 

Neutral Filters. Table 2 summarizes the results of 
many measurements made over a period of 9 months 
on the transmission values and their average devia- 
tions for the five neutral filters that are fixed in the 
instrument and moved reproducibly into position. 
The first four were always checked by the amplifier; 
whereas the very dense one was usually checked in 
comparison with another neutral filter and the ampli- 
fier. In some experiments a direct measurement of 
the dense filter transmission was made using the 
amplifier and gave good agreement with the com- 
parison measurements even though the signal to 
noise ratio was extremely low. Table 3 gives the 
transmission results under changing experimental 
viewing conditions. Also table 2 indicates that the 
1P21 and 5819 phototubes give the same transmission 
values. There is no question but that there is a slight 
nonuniformity of the transmittance over the total 
area of the neutral filter, as was seen in later experi- 
ments using a small slit to sean the filter; however, 
all ordinary measurements are carried out as above 
with fields of view of such dimensions that these small 
differences are integrated into a constant average 


set 
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value 
TABLE 2 Transmittance of neutral filte rs for different photo- 
tubes at different wavelengths of light 
1P21 phototul 5819 phototule 
Filte 
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TABLE 3. Transmittance of neutral filters at 546 


size beams and different fields of vie 


10-mm beam 3-mm beam 


Filter : , ‘ - a 
x9.9 mm 3.0X7.0 mm 5.04 mm diam 5x<9.9 mm 
| 0. 526 0. 527 0. 522 0. 552 
“ 32Y 32y 330) 323 
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2.22x<x 10-4 2 ‘19x 10-4 2? 20k 10-4 2? 24x10-4 


Table 3 indicates no differences in the transmission 
values beyond experimental error when only the 
field of view is changed. However, special attention 
should be directed to the differences which appeal 
when a beam which is smaller than the field of view 
of the photometer is used. These differing trans- 
mittances are no doubt due to the fact that a neutral 
filter with nonparallel faces or inexact perpendicular 
alinement can shift the lightbeam when placed in 
position. This shift then makes the entire lightbeam 
strike a different part of the photosensitive surface 
to give a slightly different signal. This causes no 
problem as long as the transmittance for such 
geometries is determined in the apparatus as it Is 
used. 


3.3. Optical 


There are two optical systems to consider in the 
light-scattering photometer, namely, the projection 
system and the receiving system. In this instru- 
ment both systems are telecentric, that the 
aperture stop is located on the focal plane of the 
lens so that the principal ray in the projected beam 
is at all times parallel to the optie axis... Figure 7 
shows the diagram of the system for a lens L, of 
focal length f, an aperture stop S,, and a field stop 
S,. For the projection system the dimensions are; 
f—=98.6 mm, r-—318 mm, x=48.1 mm. From these 
dimensions the angular divergence of the beam, 
represented by one-half of the angle @ in the diagram, 
and the magnification may be computed. It can 
be shown that with the dimensions in this instrument 
the angular divergence of the rays from the center 
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to the edge of the field stop differ by no more than 
2 percent in the vertical direction with the largest 
aperture stop ever used, and by much less in the 
horizontal direction. The angular divergence of the 
most divergent projected beam is calculated to be 
0.28° in the horizontal and 0.43° in the vertical 
direction, although for certain special measurements 
this beam was replaced by one having a 0.07° 
divergence in all directions. Although the diver- 
gence is difficult to measure exactly because of the 
indistinct image away from the field stop S; the 
divergence measured 250 cm away from 5S; gave 
agreement within a few percent of the calculated 
value. 

In order to determine accurately the dimensions 
of the projected beam, measurements were made in 
three ways. Visual measurements were made with 
a telescope focused on a ruled glass scale located 
on the center of rotation. Then photographs were 
taken with a camera focused on the scale, and also 
on the centering fixture used for alining the instru- 
ment. The glass scale was particularly useful 
because it could be immersed in a solvent in a light- 
scattering cell to make certain of the unchanged 
magnification when the distance x in figure 7 was 
changed slightly to accommodate the decreased 
optical distance to the center of the cell. The 
results for three beams measured at the center of 
rotation are as follows: 9.20 +0.07 by 11.97 +0.10 
mm, 3.09 +0.03 by 9.22 «0.07 mm, and 1.60 +0.03 
mm diameter. The measured magnification of 2.06 
is within 2 percent agreement of the calculated 
value. 

The receiver system may also be represented by 
figure 7 and has the dimensions x=36.50, and 
r=78 to 117 mm. The achromatic lens used in 
the receiver, however, had its principal points deter- 
mined on a nodal slide in the NBS Optics Division. 
The value for the equivalent focal length needed 
to calculate the angular acceptance of the system 
determined by 5S; and S. was 40.44 +0.10 mm for 
575-mu wavelength light. The aperture stop was 
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located on the focal plane within +0.1 mm. The 
magnification of the field stop is about 1.11. The 
entire phototube housing can be moved so that the 
optical distance r may be kept constant by changing 
the geometrical distance when large size light- 
scattering cells are used and compared to quantities 
measured in air. 

The beam uniformity has been tested in many 
wavs. When the lamp is centered on the 1- by 3- 
mm aperture stop in the projector, the uniformity 
of the 9.2-by 12.0-mm beam as measured with a 
phototube having a 0.254-mm slit in front arranged 
to accept about 2.5° semicone angle, is constant 
within +0.5 percent over 95 percent of the beam. 
In this experiment the phototube housing is advanced 
across the beam by means of a screw with a pitch of 
40 threads per inch. Evaluating more than this 
amount of the beam would require a smaller slit 
since the edge effects can be calculated to arise 
somewhere in this region. The same beam size at 
the center of rotation but having a smaller divergence 
gives the same results although the lamp is not as 
difficult to aline since less of the lamp area is required 
to fill the stop. As might be expected, the smaller 
size beams give the same uniformity. 

Alinement of the receiver and beam direction is 
necessary for these measurements, and this can only 
be achieved by the use of mechanical fixtures followed 
by optical sighting methods. The end-on phototube 
does not appear as sensitive to extremely small 
angular deviations from parallelism as the 1P21. 


The authors are indebted to P. Pfaff of the NBS 
Machine Shop for his many suggestions and excellent 
workmanship. They also-thank E. Webb of the 
NBS Optical Shop for his construction of many of 
the optical components, and his cooperation in 
fitting and mounting many of the other items. 
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Resolution of the Dissociation Constants of d,/-Malic Acid 
From O° to 50° C' 


Murray Eden ’ and Roger G. Bates 


\ method for the precise determination of the dissociation constants for the two over- 
lapping dissociation steps of a weak dibasic acid from the same set of experimental data is 
described. The procedure, based on an earlier suggestion of Speakman (J. Chem. Soc. 1940, 
855), is applied to the resolution of the thermodynamic dissociation constants of dJd-malic 
icid at intervals of five degrees from O° to 50° C. The constants were derived from measure- 
rents of the electromotive foree of cells of the tvpe 


Pt: H. g), malate buffer solution, KCILOm 0.008), AgCl: Ag. 


The 39 buffer solutions were prepared from potassium hydrogen malate and from mixtures 
of potassium hydrogen malate with perchloric acid or potassium hydroxide, and each con- 
tuined potassium chloride The first aissociation constant, A,, is 3.48 L0-* at 25° C. It 
is given as a function of absolute temperature (7') by 


: L35S8.85 - se 
lor K r 5.1382 + 0.013550 T. 
The second constant, K,, is 7.99% 10> at 25° C. It is given by 
, 1H58.53 . S-on7 
lor A , 6.2364 0.019353 7, 


The thermodynamic quantities AG®, AH’, AS°, and A( related to each dissociation step 
\ computed and compared with the corresponding constants for other acids closely 
lated structurally to malie acid 


a Introduction so successful and has been used so w idely that it may 
justifiably be termed the ‘‘conventional” emf method. 


The evaluation of the dissociation constants of a Difficulties are usually encountered, however, when 
pol basic weak acid or poly acidic base is particularly the dissociation constant of the acid or base exceeds 
difficult when the ratio of the constants for two sue- | 0.002 [33]. 
cessive dissociation steps is 100 or less. Even when The cell most commons used is a cell without 
the ratio of the first dissociation constant to the | transference, composed of a hydrogen electrode and 
second Is as high as 500, one of the constants cannot i silver-silver chloride electrode. The same cell 
be evaluated without correcting for the effect of the | was utilized in the present study. It is represented 
other, inasmuch as some simultaneous dissociation | by the following scheme 


Occurs 

In view of these difficulties, the dissociation con- Pt; H.(¢), malate solution, KC](m=0.008), 
stants of only a few dibasic or poly basic acids with Ag ‘l; Ag. (1) 
“overlapping” dissociation steps have been resolved 
In a precise way Indeed, accurate values may be The reaction taking place when current is drawn 
unobtainable by the usual experimental procedures. | from this cell is 
Consequently, little is known about the thermo- 
dynamic quantities associated with the separate over- ; H, + AgCl= Ag+ H*+Cl 


lapping lonization steps of acids of this type. 

When the hydrogen ion concentration of a buffer | Consequently, the emf / (corrected to 1 atm partial 
solution is fixed by the interaction of a single acid or | pressure of hydrogen) is related to the composition 
base with water, the dissociation constant can often | of the solution in the cell by 
be obtained with high accuracy from emf measure- 


ments by the method described by Roberts 11] hand k—E low law ( fund 9 
by Harned and Ehlers 2]. This procedure has been 2 2322627 F' OF Mtoe, OLY \JunloitH)- - 
Based lissert resented by Murray Eden in partial fulfillment In this equation, m represents molality (moles per 

1. Presented be the Division of Phys + hemistry at the 134t! kilogram of water), J ois the acuivily coeflicient on 

ga fom pe tm seg mene Heath’ Ea a nd Welles. U5} the molal scale, and £/° is the standard potential 
Bethesda 1 Ma. Institutes of ‘Health, National Heart Institute, | (standard emf) of the cell. The values for the 
3 Fiwaves tn brackets indicate the Meerature celeremecs at the end this panes latter were derived from the work of Harned and 
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Ehlers [2,4]. They have been tabulated in absolute 
volts in an earlier publication [5], together with 
values of the function 2.3026 R7T/F for temperatures 
from 0° to 60° C (absolute temperatures, 7’, from 
273.16° to 333.16° K). 

The expression relating the emf to the equilibrium 
constant, A, for the dissociation process, 


HA*"+H,O=H,O*+ A*", (3) 


where » represents the charge of the acid species 
HA, is derived by combining eq (2) with the mass- 
law expression for eq (3). The concentration of 
protons in dilute aqueous solutions is unknown, 
although it must be exceedingly small. However, 
as long as the activity of water remains substantially 
unity, the activity of protons will be proportional 
to the activity of the hydrates formed in reactions 
such as eq (3). Indeed, the hypothetical standard 
state of unit proton activity (to which £° relates) 
corresponds by definition to an actual state in which 
hydrated protons Ol hydrogen are present at 
unit activity. Hence, my in eq (2) can be identified 
With Myo” In the mass-law formulation of 
These considerations are valid only as long as the 
activity of water remains substantially unity. 
The indicated substitution gives 


1OnsS 


aD 
(>). 


log A 


log (fufeimn) 


The activity-coefficient term at finite ionic strengths 
can be estimated by various means in order to facili- 
tate the extrapolation of “apparent” values of —log 
A to infinite dilution, where the last term becomes 
zero and the true value of A is obtained. This 
method ts — elegant and completely satisfactory 
when the buffer pH is controlled by a single dis- 
sociation step of moderate — and when the 
buffer solutions are stable in contact with the hvy- 
drogen and silver-silver diestde electrodes. It ap- 
plies to polybasic acids and polyacidic bases if the 
ratios of the for successive 
steps are 10° or greater. The concentrations of the 
acid and its conjugate base are derived by correcting 


constants dissociation 


the stoichiometric molalities of these substances for 
hydrolysis and dissociation, that is, for reaction with 
the amphiprotic solvent. 

The existence of a second acidic or basic group of 


a strength comparable with that ies study com- 
plicates the determination considerably, for the 
centration term of eq (4) can no longer be computed 
simply and accurately from stoichiometri€é concen- 
trations [6, 7, 8, 9]. The quantity mg in eq (2) is 
now a function simultaneously of more than one 
dissociation equilibrium, and 


COl- 


the resolution of the 
individual constants may be a problem of consider- 
able a 

Nevertheless, it is sometimes possible to apply the 
“conventional” method suce essfully by choosing buf- 
fer ratios so that the hvdrogen ion concentration is 
controlled almost exclusively by the dissociation step 
under study a relatively small correction being made 


for the effect of the other. The procedure is an 
involved and laborious one. The dissociation con- 
stants of a number of weak dibasic acids have none- 
theless been determined through variations of this 
general procedure [6, 10, 11, 12, 13]. In this way 
the dissociation constants of o-phthalic acid (Ay/ A» 
288) have been resolved [12, 13] and the second dis- 
sociation constants of succinic acid (/t,/A,=27) |14] 
and d-tartaric acid (A,/A,=—21) |15] and the third 
dissociation constant of citric acid (A,/A,=43) |16] 
have been determined. In the last three instances a 
buffer ratio, m (neutral salt)/m (acid salt), of 
effective in keeping the hydrogen ion concentration 
low and in reducing the magnitude of the corrections 
for dissociation equilibria other than the one of 
primary concern. 

Various attempts have been made to determine the 
overlapping dissociation constants simultaneously. 
In the method suggested by Speakman [17], the first 
and second dissociation constants for a weak dibasic 
acid are derived from a single set of observations. 
The experimental data are plotted in such a way 
that a straight line is obtained. The intercept is 
Kk, Kk, and the slope determines the value of A). 
However, Speakman’s method utilizes values of 
“hydrogen ion activity” (a) from pH measurements 
and consequently does not furnish exact values of 
thermodynamic dissociation constants. 

A method for determining the product of the over- 
lapping ionization constants of a dibasic acid was 
proposed by Bates {18}. The complete resolution 
depends on an independent determination of one 
constant. The cell solutions contained known 
molalities of the acid salt of the weak dibasic acid 
and an alkali chloride. An apparent A,A, was 
plotted as a function of ionic strength and extra- 
polated to infinite dilution to eliminate activity co- 
efficient corrections and to obtain the true value of 
kK, k,. The value of this method lies in the fact that 
the molality term is insensitive to relatively large 
uncertainties jn the hydrolysis correction. In apply- 
ing this procedure to a resolution of the dissociation 
constants of tartaric acid, Bates and Canham [15] 


evaluated the second constant with the aid of 
approximate values for the first. The product 
A, Ay was then evaluated and a new, more precise, 


value of A, was computed, 


2. Method 


It seemed worthwhile to explore the possibility of 


method to permit the 


modifying the Speakman 
dissociation con- 


simultaneous determination of the 
stants of a dibasic acid with a precision comparable 
with that obtainable by the “‘conventional” method 
for a single isolated dissociation step. To this end, 
the equations were formulated in terms of fa fey 
(that is, auf<)) instead of au, since the former can be 


derived unambiguously from measurements of the 
emf of cell (1) by means of eq (2). This quantity 
will be termed P’: 

P ant, l hu i, fiw (5 
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2.1. Overlapping Dissociation 


The malate buffer solutions were all prepared from 
potassium hydrogen malate with the addition of 
perchloric acid or potassium hydroxide as desired. 
It was convenient to define two new quantities as 
follows: 
malate 


a=the amount of potassium hydrogen 


used, in moles per kg of water. 


8—the amount of strong acid used, in moles per kg | 


of water. If strong base is added instead of strong 
acid, the sign of 8 is negative. 


The potassium chloride present in all solutions at a 
molality of 0.008 must, of course, enter into the com- 
putation of the ionic strength. It does not, however, 
participate in the dissociation equilibria of malic acid, 
and it is therefore omitted in following derivation 
of the relationships among P, a, 8, Ky, and A). 

From the requirement of e ne ‘troneutralityv, we have 
the following condition: 


(K*]-+[H*+]=[ClO,-]+ [OH-] 


(HMal-]+2[Mal*], (6) 
where HMal~ and Mal® are written for the hydrogen 
malate and malate ions, respectively, and where the 
quantity in brackets represents the molality of the 
partic ular species of ion. It is convenient to define 
a quantity 2 as follows: 


B=8—[H*]+[OH-. (7) 


It will be noted that for all the malate buffer solutions 


KK | ICO, a—Pp. (S) 
From eq (6), (7), and (8), therefore, 
a— B—[HMal-]+2[Mal*]. (9) 


It can also be shown that 


a B (10) 


2|H,Mal]—-[HMal-]. 
These relationships can now be combined with the 
equations for the two thermodynamic constants, Ay 
and AY, of malic acid. In this way, an equation 
analogous to that of Speakman |17] is derived. 


2.2. Resolution of K, and K, 


Unfortunately, Speakman’s equation cannot be 
applied in an exact way to the determination of over- 
lapping dissociation constants, inasmuch as there is 
no known experimental means of measuring hydro- 
gen ion activities. The closely related function P, 
however, possesses the thermodynamic validity that 
aq lacks. With the further substitution of ag=TP fe), 
one obtains . 


| 


B Face oo 
Pp: = * = _ =P —, K,+- KK, 
ath Juvala atBb Jua-Jei 
(11) 
which may be written 
Y= Vok,+ (hk) /p (12) 
where 
— 
A ; (13) 
at+B 
and 
P2(a B) 
y , (14) 
a+B 
and in which the activity coefficient terms are 
fuyafer by 
oO (15) 
Sua 
and 
fy - > 
92-3 (16) 
Jujalc 


Apart from the correction for the concentration of 
free hydrogen or hydroxyl ions in the buffer solutions 
in the computation of B by eq (7) (see below), the 
quantities .Y and Y can be derived without difficulty 
from the stoichiometric molalities and the measured 
emf of cells of type (1). This correction is usually 
small, for the experimental conditions can often be 
so arranged that @ and 6 are large compared with 
Win and Mion. 

It can be seen that eq (12) will be linear in Y and 
) if certain conditions can be met, namely: if the 
different values of the independent variables do not 
appreciably effect the constancy of o and 1/p. 
Studies of salt effects indicate that the activity 
coefficient terms are often largely dependent on the 
ionic strength (uw) and not on the nature of the added 
salt [19, 20]. Therefore, it should be justifiable to 
consider o and p constant, within the error of the 
determination, when comparing values of .Y and Y¥ 
derived from the measurement of solutions contain- 
ing different ratios of acid salt to strong acid (or 
strong base) at constant ionic strength. This 
treatment of salt effects appears to be adequate in 
malate buffers, because a straight line was obtained 
for each value of the ionic strength when })° was 
plotted as a function of Y. The slope of each line 
is oA, and the intercept is (AV, A, /p). 

Fortunately, activity-coefficient 
same type as o are often close to unity in the regions 
of low and intermediate ionic strength (w<0.1).* At 
u-—-0.06, ¢ was found to be 0.96. Furthermore, log ¢ 
shows the expected linear variation with tonic 
strength, and consequently the curves obtained by 
plotting log (oA) against icnic strength were easily 
extrapolated to zero ionic strength. Since o@ aup- 
proaches unity at infinite dilution, the intercept for 
u—O0 is log Ay. 


terms of the 


equimolal amounts of acetic acid, sodium 
given by 0.06 uw: in acetate 
0 14 yw 21 


‘In buffer solutions composed of 
wetate, and sodium chloride, log (fuac foei/fac) i 
buffers containing potassium iodide, log ( fitac fi fae) is 
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Similarly, log A, is obtained by extrapolating log 
(KX,/p) to zero ionic strength. An activity-coefficient 
term of the charge type of p, however, cannot be 
expected to be very near unity at moderate tonic 
strengths (p was found to be about 0.6 at the highest 
ionic strengths included in this study), nor can log p 
be expected to vary linearly with ionic strength. 
In order to correct the curvature of the lines and to 
facilitate the extrapolation, estimated values of p 
(designated p° were computed from the Debve 
Hiickel formula 


log p* SAVE 7 
= > ~ , ‘ 
| 6.25 By uw 


where c.f and 2B are constants of the Debve-Hiickel 
theory |22] The quantity log p*(A./p), instead of 
log A,/p, was then plotted asa function of uw. Several 
values for the coefficient of the Byyw term in the 
denominator of eq (17) were tried, and 6.25 was 


found to furnish the best linear plot of log p*(Av,/p 
The lines obtained were readily 
strength, where p*(A,/p 


as a funetion of wp. 
extended Zero 
becomes equal to i 


to hone 


2.3. Computation of B 
In order to utilize eq (11 12 Iy 
and Ay iis described in) the preceding section, two 
other quantities, B and yu, must be known The 
former can be obtained readily from 8, the stoichi- 


or lo resolve 


ometrie molality of strong acid o1 base, by eg ri 
together with [H*] and [OH]. Inasmuch as the 
put range of malate buffers extends from about 2.5 


to 6.0, the concentration of hydroxyl ion is always 


less than | percent of my and can be neglected 
Llenece, 
3 ie] (hl | qa 
2.4. Hydrogen Ion Concentration 
Values of the molalits of hvdrogen 1On, (H |= My 
were derived directly from the emf by eq * Ff. is 


(by definition) equal to fie), the square of the mean 
ionic molal WeuIVIL coetheient of hvdrocloric aeid in 
the malate-chloride buffer solution. Values of fc, in 
pure aqueous solutions of hydrochloric acid of the 
ic strength the malate-chloride buffer 
solutions were used in the computation of [IH *] 
They were taken from Harned and Ehlers [4, 23] 


sume Lohic vs 


2.5. The Ionic Strength 


With the aid of eq (9) it can be shown that the 
pone strength of the solutions contaming nerd salt 
(KENal and potassium chloride Wi OLOUS or 


these components with added strong aeid, is given by 


u OO0NS a Niyy Mxy IS 


Inasmuch as the hydroxyl ion concentration is neg 
ligibly small 


| 


When the solutions contain potassium chloride 
potassium hydrogen malate, and potas- 


(O.0O0S m), 
sium hydroxide (stoichiometric molality of —8), the 
ionic strength ts 

u— 0.008 4+-a—28+2my+ mH.Mal 19) 


In order to calculate the ionie strength by these 
two equations, it was evaluate the 
molalities of malic acid and bivalent malate ion as 
well as that of hydrogen ion. These quantities can 
only be obtained through successive approximations, 
as values of the ionic strength and also of the ap- 
proximate dissociation constants must be used in 
their computation. However, all were small rela- 
tive to a, and usually to 8 also, so that the approxi- 
mations converged rapidly 


necessary to 


The molality of malate ion, needed to calculate the 


ionic strength by eq (18), can be expressed by 
AK - 
Mii Mal »» 3 Min Mii Mal 2U) 
ol - 
with the aid of eq (5), (Va), and (16 In the solu- 


tions of acid salt, to which no strong acid or strong 
base has been added, OB 1s zero and consequently 


Ky kK. Min 


5)? KK " 


Mints 


In the solutions to which perchloric acid had been 


added, 


Oo was large compared with Win and Wingai™; 


STs] 
, pl’ 


Mints 21a 


When the solutions contained added strong alkali, 


it Was necessary to estimate the molality of free 

malic acid in order to compute the tonic strength by 
eq (19 Krom eq (5), (7a 10), and (15 

|P ok, a re} Min 90 

ee I+2PicK,)) 22 


alkali had been added to the solution 


When strong 


of acid salt, 2/7? ¢A, was much less than unity. 
Furthermore, o is nearly unity at low and inter- 
mediate tonic strengths: henee, 
. pP 
Min > cr 3 iy 22a) 
a 
The accurate evaluation of mig oar, aNd Maggy? as 


described above Is quite evidenthy it laborious pro- 
cedure. However, neither of these quantities makes 
a large contribution to the tonic strength obtained 
yy (18) and (19 and therefore accurate values 
are unnecessary. Consequently, the following pro- 
cedure was followed: A> preliminary value of the 
lOnLLC strength was derived with the omission of the 
last two terms of eq (18) and (19), and a first approxi- 
mation to m_ was made, which led to a more accurate 


eq 
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value of uw by eq (18) and (19), the last term of each 
of these equations being omitted. 

The first approximation to mya and My,Mai Was 
then made by eq (21), (21a), or (22% the choice 
being dependent on the composition of the buffer 
solution. For this purpose, the values of kK, and K, 
could be obtained with sufficient accuracy from the 
literature, from titration curves, or from a_pre- 
liminary plot of Y as a function of VY. The ionic 
strength was now recalculated, a new value of my 
obtained, and the process repeated until a consistent 
set of values resulted 


3. Experimental Procedures 


Potassium hydrogen malate was prepared from 
reagent-grade dJ-malic acid and potassium bi- 
carbonate A solution of the salt was decolorized 
at 60° © with activated charcoal. The filtrate was 
cooled in a water-ice bath and the crude salt precipt- 
tated by the addition of ethanol. This precipitate 
was recrystallized five times from ethanol-water 
mixtures. The salt collected after the fifth crystal- 
lization was dried overnight at 130° C in air. 

The potassium hydrogen malate was assayed by 
titration with the standard solution of potassium 
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hydroxide described below. After the first recry stal- 
lization, the dried salt was found to assay 99.69 
percent potassium hydrogen malate. After the 
third recrystallization, the assay was 99.99 percent. 

Potassium chloride was prepared according to the 
procedure outlined by Pinching and Bates [24]. 

The standard carbonate-free solution of potassium 
hvdroxide (0.1 molal) was standardized against 
potassium hydrogen phthalate (NBS Standard 
Sample 84b) in a carbon dioxide-free atmosphere. 
Weight burettes were used. This alkali solution was 
used both for titrations and for the preparation of 
the cell solutions. 

Reagent-grade perchloric acid, 70 percent by 
weight, was diluted with conductivity water to make 
a solution of molality about 0.1. The acid solution 
was standardized against the potassium hydroxide 
solution. It was tested for the presence of chloride 
with silver nitrate solution according to the pro- 
cedure described by Rosin [25] and was found to 
contain less than 0.001 percent. 

Kight stock solutions of malate buffers were pre- 
pared. Additional cell solutions were prepared by 
diluting these with an 0.008-m solution of potassium 
chloride. The compositions of the solutions appear 
in the second and third columns of table 1. 


rived from emf measurements of cell (1 
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The apperatus and the experimental procedures | 4 ..o¢%———-—— a a u 
emploved in this work were essentially as deseribed a P 
by Bates and Acree [26]. The temperature of the ° “ 
bath containing the cells was maintained within s 

-().02° C of the nominal temperatures. The temper- Ta 2° 
ature sequence used, with a few exceptions, was 25°, | Pd e 
15°. 5°. 0°. 10°. 20°, 30°. 40°, 50°, 45°, 35°, and 25 a 4 5 e 
(.. Individual emf readings were discarded if the Pa ond 
cell potential varied by more than 0.2 mv from the y, as = 
smooth curve obtained by plotting emf against 4 " ae e ~ 
temperature. All measured values for a given solu- 
tion were discarded if the final value of 25° C differed ™ ° 
from the initial value by more than 0.3 my. o 

The cells usually attained their equilibrium values ; . 
in about 30 min after each change in temperature. ° 
Equilibrium was judged to have been reached when e 
readings of the same set of electrodes taken 10 min ° | 
apart differed by 0.03 mv or less. 

4. Results and Calculations 

The values of —log (fafevra) for each solution 
from 0° to 50° C are assembled in table 1. The | oh gy : 
solutions were composed of potassium hydrogen ¢/ 
malate (molality a) and potassium chloride (molality é 
0.008), and some of them also contained perchloric 
acid or potassium hydroxide (molality +B). <A 
positive value of 8 signifie ‘s the presence of perchloric | Frcure 1. Values of X from 0° to 50° C plotted as a function 
acid, whereas a negative value indicates that potas- ilies of tonic strength 
sium hydroxide had been added. Each figure for 

log (fafeyma) represents the mean value derived 
from the two pairs of electrodes in the same cell. 

“ 


4.1. Calculation of K;, 











The quantities VY and ), computed by means of |... -« ;, , a 
eq (13) and (14), were plotted as a function of ionic . $° 
strength. Different curves were obtained for each 3° 
series of solutions and for the same series at different s a 
temperatures. Two representative plots are shown - e: 
in figures 1 and 2 The values of Y and Y corre- e ov 
sponding to tone strengths of 0.06, 0.05. 0.04. 0.08. * e 
0.025, 0.02, and 0.014 were obtained by interpolation . od 
on each curve { va . 

A second series of graphs was prepared, in which * ° ° 
Y was plotted as a function of .V at each of the seven : 
values of the ionic strength and at each of the 11 PY : . 
temperatures. These plots consisted of straight W///, 
lines, as shown in figure 3. The position of each line 3 
was determined by the method of least squares. /g P 
The logarithm of the slope (@A,) of these lines was ; i 
next plotted as a function of ionic strength at each 4 
temperature (compare fig. 4). The sorte res , bs 
straight line drawn through the points was extended 
to w=0, where log ck,=log A,. The mean devia- 
tion of the individual values of log cA, from the line 
was +0.0010 unit. The values of log A, obtained Pi 
in this way are summarized in table 2. Z 

4.2. Calculation of K, 

The second dissociation constant was derived from 
the intercepts (A, A p) of the plots of ) with respect Ficgurk 2. Values of Y from 0° to 50° ¢ plotted as a function 
to \. If the experimental errors for all of the meas- Series 1 5 a nailenanaa cal 

166 
- 

















ured points were about equal in magnitude, it is to TABLE 2. Summary of log Ky, and kK, from 0° to 50° ¢ 
be expected that the best values of the intercept 
would be obtained by the shortest extrapolations. 


: : - ; ~ remperature —log K hk, X10 
For this reason, only the data derived from the four 

series (4 through 7) composed of the least acidic C 

solutions (8<0) with the smallest values of VY and ° . = oe 

)” were used to derive Ky. The value of kK, K, p cor- 10 3. 494 3. 20 

° » . . 15 3. 482 2v 

responding to each of these experimental points was 2) 5 ar? 37 

25 3. 458, 3. 48 

+0) 4, 452, $53 

35 4, 446) 4. OF 

7 1) $ 444 3. 00 

4 3. 446 3 5S 

v1 $. 445» 3. 5Y 

















e . 
+ os e 
; ee® 
———. em “Ee a wa ve; 
a oC ef e* ” - . 
i> 5 a ° = 
> , e e* e@«e 
pence . a re er 
. « 550 7. « 
q 2 +. e s a 
- + e 
5.16 ~ e@ @e 
} 2. ____ gy _________ - @ ____ ——-~9-o- - 
$14 59° @ e ee? °. 
e a 
04 ¢ 


Ficure 5. Plots of —log p*(Ko/p) at 0°, 25°, and 50° Cas a 
function of tonic strength. 


computed from the individual values of .Y and Y, 
together with the slope of the appropriate .Y, )° plot 
determined by the method of least squares. 





The values of log p*(A2/p) were independent of 

/ ionic strength, within experimental error, as can be 
seen in figure 5. It appears, therefore, that p* caleu- 

i + | lated by eq (17) is a very good representation of p 
' ' | at these ionic strengths, and it follows that p*(Ay/p) 
is equal to Ay. The mean values of log Ay at each 











( 3 lots of Y sus 6° C. . . . “1 ; ” hy} 
FiGuR} Plots of Y versus X at 25° C. temperature are summarized in table 3. An arbi- 
From left to right, the lines correspond to ionic strengths of 0.06, 0.05, 0.04, 0.08, | trary weight of 2 was assigned to those values com, 
125, 0.020, and 0.014 I'he axis of ordinates has been transposed to the right by . ‘ m 


for each successive member | puted from the 10 solutions (series 4) which contained 
neither strong acid nor strong base. For these solu- 
tions, 8 was 0. The values of Y were consequently 


‘to small, and errors in establishing the slope of the , 
ee Y line had only a slight influence on the intercepts 
bd ie ee btained 
a oe obtamed., 
a « — es j 





TABLI 3 Summary of log kK, and k, from 0° ta 50° C 
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4.3. Dissociation Constants and Temperature 


The values of log A, and log A, can be expressed by 
the equation suggested by Harned and Robinson 27). 


namely 


—log K 


1. BR. and C’ were determined by the 


The parameters . 
and the following results 


method of least 
were obtained: 


squares, 


1358.85 


log K, , -5.1382 +0.013550 7’ 24 


1658.53 


T 


log K. 6.2364+-0.019353 7, 


where 7’ is the temperature on the Kelvin seale. 
The average differences between the “observed”’ 
values of log A, and log A, at the 11 temperatures 
and the corresponding values calculated by eq (24) 
and (25) were 0.0013 unit and 0.0011 unit, 


respect i\ ely ; 


4.4. Comparison With Earlier Work 


A considerable number of determinations of the 
dissociation constants of malic acid has been made 
in the last half century. The values obtained by 
these earlier workers have been assembled in table 4. 
The columns headed Ay and kK, list values of the 
“concentration constants’, uncorrected for activity 


coefficients. The columns headed A, and A, give 


Summary of reported 


values for the thermodynamic constants. Unless 
the author made a definite statement that activity 
coefficients had been used to correct the equilibrium 
equations, it was assumed that the dissociation con- 
stants had been calculated on a concentration basis. 
The values of the “concentration” dissociation con- 
stants AY and Ky, listed in table 4 are uniformly 
higher than the thermodynamic dissociation con- 
stants. This is to be expected, since the activity- 
coefficient correction would decrease the value of the 
dissociation constant in dilute solutions. 

It is clear from table 4 that agreement 
the results of various workers is much better for the 
first constant than for the constant. All 
of the earlier values of A, are somewhat lower than 
those found in this investigation. Only the electro- 
metric determination of A, by Larsson [28] agrees 
well with the value reported here. No statement 
can be made about the relative accuracy of these 
earlier determinations of Ay. 


between 


second 


5. Thermodynamics of the Dissociation 
Processes 


The thermodynamic quantities associated with 
the two dissociation steps were computed by the 
following equations: 

AG 2.3026 RCA} BT. Cr 
AH? =2.3026 R(A—CT?*) 

AS 2.3026 R(B+-2CT) 
AC’ 2.3026 R(2CT) 


26) 
27) 
28) 
29) 
where FR is the gas constant, 8.31439 } deg~'mole 

The thermodynamic constants are summarized in 


table 5. 
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TABLE 5 Thermodynamit quantities for the 


malic acid and hydrogen malate ion from O° to 50° C 


Temperature ; AH AS 


HoMal*+H*++H Mal 


Process 


mole j deq=' mole deg™' mole 
6, 662 43.4 
5, 946 15.0 
2149 4s. 6) 
477 o1.2 
«23 53. 7 
, eh 3 
4 17 6 SU 
61.5 
4 


210 
7. 650 

, 120 

590 
Zz two 
24. 610 
30, 120 
30, 680 
31, 250 
31, 840 


6. Discussion 


It is of interest to compare the dissociation con- 
stants and associated thermodynamic quantities for 
the three structurally related acids, succinic, malic 
and tartaric. The constants A, B, and C of eq (23 
for the dissociation constants of these three acids 
are as follows: 


A B ( 
Sueccinie - id 1206. 25 3. 3266 0. OLL697 
Malice acid 1358. 85 5. 1382 013550 
Tartarie acid 1525. 59 6. 6558 015336 


1679 1 3 5. 7043 
1658. 53 6. 2364 
L765. 35 7. 8015 


. OL91L53 
O19353 
019276 


Suecinie acid 
Malice acid 
Tartaric acid 


thermodynamic data derived 


A comparison of the 
from these constants is presented in table 6. 
It is noted first of all that progressive substitution 


of hydroxy groups on the two middle carbons 
enhances the strength of both acidic groups. This 
effect is the expected one. Thus, the value of log 


TABLE 6 Compa 


*Ssuccini 
COOH-CHeCH 

Malic 
COOH-CHOH-CH 
lartari 
COOH-CHOH-CHOH-COOH 


COOH 


CooH 


. Pine hi 
b Bates an 


495640 


dissociation of 





son of the the / modynami constants for the two dissoc ration ste ps Of Succinie 


K for glycolic acid [43] is 0.9 unit lower than that 
for acetic acid [2, 44], and that for lactic acid [45] 
is lower by 1.0 unit than that for propionic acid 
146). The decreases are somewhat smaller for both 
steps in the dissociation of the dibasic acids. 

The interpretation of the changes of entropy and 
heat capacity associated with the dissociation of 
weak dibasic acids is not simple. Statistical, electro- 
static, and structural factors have been discussed by 
Everett and his coworkers [47, 48]. It is likely that 
orientation of solvent molecules in the electrostatic 
field adjacent to the ions plays an important role in 
determining the sign and to some extent the magni- 
tude of the changes of entropy and heat capacity 
when an acid dissociates. Structural features lead- 
ing to “‘sereening”’ or solvent exclusion have an oppo- 
site effect. When the dissociation produces an in- 
crease of charge, as it does with uncharged acids and 
with acid anions, the increased solvent orientation 
and attendant immobilization of solvent molecules 
should bring about a decrease in both heat capacity 
and entropy. Indeed, and ACS for acids of 
these types are almost invariably negative. 

Although solvent orientation is recognized as a 
factor of primary importance, other factors often out- 
weigh its effect, as shown by the contrary variation 
of AS° and AC with chain length in a series of 
acids of the same charge type [47], and by the fact 
that the values of AS®° for the dissociation of some 
cationic acids are negative while ACS is positive 
(51, 52]. At the present time, even qualitative 
predic tions are often of little value. 

It is nevertheless worthwhile to compare the 
observed values of these thermodynamic quantities 
for acids of like structure, to note the similarities, and 
to attempt toe xpl iin the differences. The-constane y 
of AC? for the second dissociation ste ps is perhaps 
the most striking feature of the comparison given in 
table 6. The regular increase of AS® from succinic 
acid to tartaric acid suggests that progressive sub- 
stitution of the two carbon atoms lying between the 
carboxyl groups is accompanied by an increased ex- 
clusion of solvent molecules. The concomitant de- 
crease of ACS for the first dissociation step is also 
regular, but the reason for a change in this direction 
is not readily apparent. 


See, for example, references [49] and [50 


, malic, and tartaric 


second sté p 
AH 


mole ‘ mole deg: mole 
0 109 
1181 101 


a4) SU) 
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Capacity Requirement of a Mail Sorting Device’ 
B. K. Bender and A. J. Goldman 


A mathematieal model of a sorting device 


suggested by S. Henig is considered. The 


relevant parameters are r (the number of destinations for the mail) and k (the number of 
letters entering the device during each cycle 


will never jam if its capacity is at least (7 
capacity requirement, 
viously known estimate 77k. 


1. Introduction 


We deal with a hight, idealized mathematical 
model of a mail sorting device suggested by S. Henig 
(NBS Electronic Instrumentation Section). For our 
purposes, | he operat ion of the device can be described 2 
as follows. Mail (to any of r¢ destinations) enters 
the system. After & letters have entered, the device 
‘asks itself’ to which destination it contains the 
most letters. All letters to this predominant destina- 
tion are then dropped out of the device, another & 
letters enter, and the process continues, 

The device is said to jam if, after a dropout, it is 
still so “full” that entrance of the next & letters would 
cause an ‘overflow’. We will determine the mini- 
mum Capa itv require “l of the device to ensure that 
no jamming oceurs. This is of course equivalent to 
determining (in terms of r and &) the marsimum pos- 
sible contents of the device under the dropout rule 
described above, and we shall work with this alter- 
native version of the problem. The derivation in- 
volves nothing more complicated than counting up 
letters, so that nonmathematician readers should be 
able to follow the arguments. 


2. Statement of Results 


In order to describe our results, it is convenient 


to define 


r(t) number of letters in the device just before the 
fth dropout. (1 


has shown that 


A. Bruce Clarke ' 


r(t) <?°k, for all ¢, (2 

which proves that the number of letters in the device 
never exceeds 7°k, so that the device will never jam 
if its capacity is 7k or more. Clarke (see footnote 4) 
remarks that “this upper bound for x(f) is clearly much 
too crude to be of any practical use;”’ for example, for 


Part of a project supported by the Post Office Department, Office of Research 


nd Engineerin 


2 The physical devi “dropout rules” other than the 


1 also Operate under 


one described below 


4 rule for breakir **ties’’ between destinations is also required 
‘A. Bruce Clarke, A mathematical model] for the Henig sorting machine, partly 
ibstracted in Ann. Math, Statistics 29, 622 (195s 
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of operation). It is shown that the device 
\(k—1)-4 
and for realistic values of r and & is significantly 


k letters; this is the lowest possible 


less than the pre- 


the Philadelphia-mail-type data (with r=200, k=3) 
used for numerical illustration, it only informs us 
that the device will never jam if its capacity Is 
120,000 or more. 

As was pointed out by L. 8S. Joel (NBS Computa- 
tion Laboratory), the estimate (2) can be ve rv much 
reduced using rather simple arguments. It will in 
fact be shown that 

r(t)<(r—1)(k—1)+4&, for all ¢, (3) 
which proves that the number of letters in the device 
can never exceed (r—1)(k—1)+4k, so that a capacity 
of at least (r—1)(k—1)+-- will ensure that the device 
never jams. For the Philadelphia-tvpe data, this 
shows that a capacity of 401 (rather than 120,000) is 
sufficient to prevent jamming. The proof of inequal- 
itv (3) will be given in section 3. 

In addition, it will be shown that the device can 
jam if its capacity is /ess than (r—1)(k—1) +h; that 
is, it is possible (under the dropout rule described 
above) for the number of letters in the device actually 
to reach (r—1)(k—1)+k. In other words, the esti- 
mate in (3) is not only an upper bound for «(f) but is 
the maximum possible value for «(f): 


r(t) (r—1)(kK—1) +k (4) 


max poss 


This will also be proved in section 3. Section 4 con- 
tains some comments on the probabilistic aspects of 
the problem. 


3. Proofs 


3.1. Proof of Inequality (3) 


It ‘s assumed, of course, that the device had fewer 
than ( 1)(A—1)+-4k letters in it at the start of opera- 
tion. An “indirect proof’ will be used; that is, we 
can tentatively suppose that (3) is false and show 
that this supposition is untenable. 
typ of t (ie 
Since (3) is false 
we have 


If (3) is false, there must be a first value 

a first dropout) for which it fails. 
for tt) but true for t=t—1. 
r(t,) >(r—-1) (k—1) +f, (5 


r(tg—1) << (r—1)(kK—1) +k. (6) 








We continue the argument by pointing out two 
facts about the contents of the device directly after the 
(t¢ —1 )-st dropout: 

(a) There are more than (7 
device, and 

(b) There are at 
some one destination. 

The truth of (a) follows directly from (5) and the 
fact that only & letters were added between the 
(fg—1)-st and f-th dropouts. All letters to some 
destination left the device in the (f,—1)-st dropout, 
so that if (b) were false the device would contain at 
1) letters to each of at most (7—1) destina- 
tions, at most (r—1)(k—1) letters in all, con- 
tradicting (a). Thus (b) is true. 

From (b) and the dropout rule it can be seen that 
the (4 | )-st dropout consisted of at least & letters, 
so, examining the contents of the device just before 


-1)(k—1) letters in the 


least & letters in the device to 


most (/ 
or 


the (f6—1)-st dropout, we have 
(number of letters to predominant dest > he 
and by (a), 
(number of letters to other dest.’s) >(r—1) )k—1); 
the last two inequalities together vield 
r(tys—1) >(r— 1) (k—1) +f, 7) 


which contradicts (6). Thus the supposition that (3 
is false Is untenable, so (3) must be true 


3.2. Proof of Equation (4) 


We already know, by (3), that the device can never 
contain more than (7 | \(k 1)- k letters Thus, in 
order to prove | 1), it suffices to exhibit a sequence of 
possible events leading to a situation in which the 
device contains exactly (r—1)(kK—1)+hk There are 
any such sequences, and the one described below is 
not necessarily the shortest one. 

The sequence will be constructed in two stages, 
beginning with the device empty at f=0. At the 
end of stage 1, the device will contain /—1 letters to 
each of 2 destinations and £—2 letters to each of the 
other r—2 destinations. At the end of stage 2, the 
device will contain exactly (7—1)(k—1)+* letters 
and the proof of equation (4) will have 
completed 

Stage 1: We suppose that the ~ destinations are 
numbered Ith some Way. The sequence of possible 
events begins as follows: ° The first & letters entering 
the device consist of one letter to the first destina- 
tion and & 1 letters to the rth destination; these 
last (A—1) letters then leave the device in the first 
dropout. The second & letters to enter the device 
consist of one letter to the second destination and 


been 


k}—1 letters to the rth destination, ete. After the 
(7 1)-st dropout, the device contains no letters to 
the rth destination and one letter to each other 


destination. The rth set of & letters to enter the 
device consists of one letter to the first destination 


For a more forma leseriptior if the f wing proce et Rit 


_ for t=(r—1 


and k—1 letters to the rth destination, and the 
process continues as before; after the (r—1)(k—2)-nd 
dropout the device contains no letters to the rth 
destination and k—2 letters to each of the (r—1 
other destinations. 

The & letters entering the device between the 
(r—1)(kK—2)-nd and ((r—1)(k—2)+1)-st dropouts 
consist of k—1 letters to the rth destination and one 
letter to some other destination; the device now 
contains kK—1 letters to each of 2 destinations and 
k-—2 letters to each of the other r—2 destinations. 
In other words, if we define 


number of destinations to which there are) 
kk}—1 letters in the device just before the? (S) 
fth dropout, ) 


nit) 


then we have 


n((r—1)(k—2)+1)=2. (9) 
The ((r—1)(k—2)+1)-st dropout will consist of 
k—1 letters to whichever of the 2. destinations 


mentioned above is “preferred” by whatever tie- 
breaking rule is emploved. After this dropout the 
device will contain no letters to some one destina- 
tion, £—1 letters to some other destination, and k—2 
letters to each of the r—2 remaining destinations. 

Stage ?: In this Stage matters will be arranged so 
that 

n(t) nit 1)+] (10 

(k—-2)+2<t<(r—1)(k—-1 


for ' (r—] 


i.e., so that after each dropout there are k | letters 
in the device to one destination than before 
If this is done, then by (9) and repeated application 
of (10) we have 


Miore 


n((r—1) (k—1)) =n ((r—1) (k#—2) + (r—1)) =, 

so that just before the {7—1)(k—1)-st dropout thi 
device contains & | letters to each of the 7 destina- 
tions. Then ’&—1 of these r(/— 1) letters are dropped 
leaving (r—1)(k—1) in all, and & new letters enter; 
the device now contains (7 1)(k—1)-+-# letters, and 
equation (4) is proved. 

To achieve (10) we choose the & letters entering 
the device between the (¢ l)-st and ?fth dropouts, 
(k—2)+2, (r—1)(k-—2)4+3, . a 
(7 Ly(hk—1 to consist of k—1 letters to the desti- 
nation just previously dropped out, and one letter 
to a destination which previously had /—2 letters 
to it in the device. 


4. Probabilistic Aspects of the Capacity 
Problem 


We have been concerned above with choosing the 
capacity of the device so that jamming never occurs 
In practice we might well have the more modest 
aim of keeping the frequency of jamming down to 
some tolerable level, and so could ret by with il 
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smaller capacity than that indicated by eq (4). On 
the one hand this leads to the “policy”? question of 
what jamming frequencies are to be considered 
intolerable; on the other hand, it leads to the purely 
mathematical problem of what capacity is needed in 
order to keep jamming below these frequencies. 
This last problem seems very difficult. As was 
pointed out (unpublished memorandum, July 1956) 
R. Rosenblatt (NBS Statistical Engineering 
Laboratory), the behavior of the device constitutes 
a stochastic process (indeed, a Markov chain) 
governed by the probabilistic distribution of mail by 
destinations. The possible number of “states” of 
the device is very large (for realistic values of 7 and 
k-), and the creation of analytic techniques for the 
treatment of Markov chains with a great many 
states appears to be one of the “underdeveloped 
areas’’ of applied mathematics. We know that the 
specific Markov chain in the current problem has a 
probabilistic “‘steady state’? toward which the device 
“settles down,” and the average number E(x) of 
letters in the device before dropouts, in this steady 


by J. 


state, might serve as a very rough initial guide ® to 
the capacity required of the device. Even the de- 
termination of #(2) seems far from easy; we may 
note that the upper bound derived for A(x) by a 
complicated probabilistic argument reference 
cited in footnote 4) yields the result A(2)<545 for 
the particular illustrative data used, whereas our 
eq (4) (which was not aimed at estimating /(r) and 
does not use the probabilistic distribution of mail by 
‘an be combined with the inequality 


{sec 


destinations) 
k(x) sri imax poss 


to \ ield the sharper result (2) <401 for the same 
data. 


unplitude oscillations around this average must also 


Che frequency of lara: 
required ot the 


be considered; at the very least, some information is 
deviation of the number of letters in the device 


Wasuineton, January 5, 1959. 
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Effect of Internal Radiant Heat Transfer on Temperature 
Distribution, Thermal Stress, and Deflection in Box Beams - 


Stanley Gcodman, Stanton B. Russell, and Charles E. Noble 


Thirteen box beams were analyzed to determine the effeet of internal radiant 
transfer on temperature distribution, thermal stress, and deflection. 
check of temperature distribution and beam deflection was made for one case. 
F, 
of the beam had an appreciable effect on the cover (but not the web 
uneven greater extent, on the beam deflection. 


mum beam temperatures above 700° to 900 


heat 
\ rough experimental 
For mavxi- 
change in emittance of the interior surfaces 
temperatures and to 
At maximum beam temperatures of 1,200° F, 


an inerease in inter‘or surface emittance appreciably relieved the maximum thermal stress 


1. Introduction 


The subjection of aireraft and other structures to 
increasingly high temperatures gives rise to a need 
for predicting the transient temperature distribution 
in the structures, and the resulting thermal stresses 
and deformations. This knowledge will make pos- 
sible the rational design of the structure, and the 
devising of methods of alleviating thermal stresses 
and deflections. 

As the temperature of a structure increases, radia- 
tion becomes more important as a mode of heat 
transfer [1]. The influence of radiant heat transfer 


on the temperature distribution of the structure is | 


affected by structural geometry, heating rate, and 
thermal properties of the beam material. The 
purpose of this investigation was to determine the 
effect of changing the emittance of the interior sur- 
faces of a box beam on temperature distribution, 
maximum thermal and beam deflection in 
heated box beams. 


stress, 


2. Structures and Material 


Thirteen analyses were made of modifications of 
three basic types of box beams. The beam cross 
sections are shown in the figure at the top of table 1; 
wall thicknesses are indicated by the ratios listed in 
columns 2 and 3 of the table. One type (cases 1, 2, 
and 3) is a thick-walled beam, similar to a beam used 
in the experimental work. The other types are typi- 
cal of beams used in aircraft structures; one type is 
relatively thick-walled (cases 4 to 8) and the other 
relatively thin-walled (cases 9 to 13). 

Thermal and elastic properties were taken as 
those of the type 302 stainless-steel test beam. Ther- 
mal conductivity and specific heat were approxi- 
mated by the linear relationships: 


k= 7.08 +-0.0043 7 (1) 

c= 0.106+ 0.0000257 7, (2) 

rt work w lucted at the National Bureau of Standards under the 
ponsorsh ip and with the financial assistance of tre Office of Naval Research 

Figure n kets indicate the literature references at the end of this paper 


where k is thermal conductivity (Btu/hr ft ° F); 
cis specific heat (Btu/Ib ° F); and 7 is absolute tem- 
perature (° R). Hemispherical total emittance was 
taken as 0, 0.35 (emittance of the test beam, cases 
1 to 3) and 1, as indicated in columns 4 and 5 of 
table 1. 

Young’s modulus and the coefficient 
thermal expansion were represented by 


of linear 


k—=31.0™ 10°— 0.0064 7 ™ 10° 


a=8.29 10-*+0.00137 7 107°, (4) 


where / is Young’s modulus (lb/in.*), @ is coefficient 
of linear thermal expansion (in./in. ° F), and 7 
absolute temperature (° R). 


is 


3. Boundary Conditions 


The beams were assumed to be subjected to a 
uniform heat input along one cover in a rarefied 
atmosphere. The heating rates are listed in columns 
6 and 7 of table 1.) The initial beam temperature 
was assumed to be 78° F. For the test) beam 
(cases 1, 2, 3), the boundary condition was taken 
as the heated cover temperature (figs. 2 and 7) 
observed in an experiment. For all other cases, the 
boundary condition was taken as the assumed heat 
input to the cover. All elements of the heated 
cover were assumed to remain constant at 1,200° F 
after reaching that temperature. 


4. Method of Analysis 
4.1. Calculation of Temperature Distribution 


Analysis was made by a numerical method similar 
to that deseribed by Dusinberre [2]. The beam 
cross section was divided into analysis elements, 
and a heat-balance equation was set up for each 
analysis element and solved for the element tem- 
perature after a short time interval, Aé@ The 
process was repeated for successive short-time 
intervals using the new element temperatures as 
the starting point of the next computation. It 
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TABLE 1 ( omp led cov lem pe ratures, maximum thermal stress and marimum deflection in box beam heated along one cove 
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was assumed that) surface emittance remained | the external heat transfer was svmmetrical wit! 


constant and that beam surfaces reflected diffusely respect to the left and right halves of the beam 
and were gray, Le., total emittance and absorptance | cross section, the right half of the beam was repre 
are equal at all beam Lemperatures Craseous heat sented by i perfect but diffuse reflector. The 
transfer was neglected reflector was divided into four regions The element 


One-half of the cross section of the symmetrical | configuration for cases 1, 2, and 3, shown in figure 1 
beam was divided into 12 analysis elements since Is typical of all cases 
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The heat balance equation is as follows: 


anne _ Aé _ a 
T..= T,.+-=+—*— (kem—1.m) (Hm) (Tn ~1 — Tn) 
a 4a 


(Kem. m1) (Fu) (Tna-s— Tn) + €n6Am(T4—TS) 
Al [F'n (T{—T3) + Fn.o(T1—T3) 


Frnio(T4—TA)|+Wr >> 


where 


T,, is the temperature (°R) at the center of element 
m after time interval Aé@ 

V,, is the volume (ft*®) of element m 

7,, is the temperature (° R) at the center of 
element m at the beginning of the time interval 
A@ 

p is the density of the material (Ib/ft*) 

(’,, is the specific heat of element m (Btu/lb ° F), 
represented by a linear function of 7,7. 

ere * are the thermal conductivities of 
the beam metal at the mean temperatures of the 
elements m and m+1, and m and m—1, 
respectively (Btu/hr ft ° F) 

HI,,, J, are ratios of element contact 
distances between element centers for element 

respec- 


areas to 


m and adjacent elements m—1, m+-1, 
tively (ft) 

e, is hemispherical 
surface m 


total emittance of element 


o is the Stefan-Boltzmann constant (0.1713 * 107°) 
in Btu ft? hr ° Rt 

A,, is exterior surface area ({t°) of element m 

T, is ambient temperature (° R) 

Aj, is interior surface area (ft®) of element m 

Ff,» is an over-all radiant heat interchange factor 


for net radiant heat exchange between a gray 
surface of element m and a gray 
element 7 It includes the effect of direct and 
all reflected radiation. 

He, is the external heat input (Btu/hr). 


The first two terms inside the brackets on the 
right-hand side of eq (5) represent conductive heat 
transfer between element m and adjacent elements; 
the third term represents radiant heat transfer with 
the outside environment, the fourth term represents 
internal radiant heat exchange between element m 
and the rest of the beam interior; the fifth term 
represents the assumed heat input boundary condi- 
tion. For cases 1, 2, and 3 in which the boundary 
condition was taken as the observed heated cover 
temperature history of the test beam, the temper- 
ature of each heated cover element was expressed by 
four linear functions of time covering successive 
time intervals. 

To compute the gravbody radiant heat exchange 
factors, F, blackbody radiant heat exchange factors 
were computed for all two-element combinations of 
the 16 internal surfaces of the analysis elements and 
(fictitious) reflector using the methods of reference 
[3], and radiation geometry. Values of the 256 F's 


surface of 


were then obtained from the matrix relationships of 
reference [4]. The computations were performed on 
SEAC using an existing code for inverting the 
matrices. 

The temperature-distribution computations were 
computed on SEAC. Temperatures at the mid- 
points of the 12 analysis elements of the beam were 
printed out at regular time intervals which were 
whole number multiples of Adé. 

In order to minimize computing machine running 
time, it was desirable to use the largest value of the 
time interval, A@, consistent with adequate accuracy 
of solution. The time interval must be sufficiently 
short that during the interval (1) The initial element 
temperatures can be used with negligible error and 
(2) for conductive heat transfer; the effect of any 
nonadjacent element on a given element is negligible. 
Satisfactory values of A@ were obtained by trial. A 
portion of the temperature distribution history was 
computed several times using successively larger 
values of A@. The largest value of A@ which gave a 
temperature distribution history negligibly different 
from that obtained with the smallest value of A@ was 
used to compute the complete temperature distribu- 
tion history. Values of A@ found satisfactory by this 
method ranged from 0.72 see for rapidly heated 
beams to 1.8 see for slowly heated beams. 

Some error was generated in the beam temperature 
history computations by representing one-half of the 
beam cross section by 2 perfect, diffuse rather than 
specular, reflector. As a check on the magnitude of 
this error, radiant heat transfer rates were computed 
for a simple symmetrical case using first the entire 
beam cross section in the computations, and then an 
equivalent beam consisting of one-half of the beam 
cross section and a diffuse reflector. An infinitely 
long box beam of rectangular cross section divided 
into six analysis elements was used. The beam was 
2-in. deep, 8-in. wide, and emissivity was 0.35. The 
vertical (2-in.) walls of the beam were at temper- 
ature 7) and the horizontal walls at absolute zero. 
The rate of radiant heat transfer to the left half of 
the upper horizontal wall was 


qg—0.1354e7} for original beam 


q O.138907) for equivalent beam 
with diffuse reflector, (6) 


where gisrate of radiant heat transfer (Btu hr ft*); 
o is the Stefan-Boltzmann constant (0.1713 107° 
Btu (ft? hr ° R*)); and 7) is vertical wall tempera- 
ture (° R). It was concluded that the error due to 
use of a diffuse rather than a specular reflector was 
small 


4.2. Computation of Thermal Stress and Deflection 


Thermal stresses and deflections in the beam were 
computed by the methods described in reference [5]. 
Integrations over the beam section were 
performed numerically using finite elements identical 


cCTOss 
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with those used in the temperature distribution 
computations (fig. 1). Possible effects of vielding, 


creep, and buckling were neglected. 


5. Test of Beam 


The test specimen, corresponding to case 2, table 
1, was a box beam constructed of 0.128-in.-thick 
type 502 stainless-steel sheet (fig. 1) Two pieces 
of sheet were each bent into identical channel shapes 
and butt welded together lengthwise to form the 
2 in. by 5 in. by 30 in. beam 

The beam was heated in a vacuum chamber whose 
pressure maintained at 4 mm of mercury, a 
pressure low enough to eliminate appreciable con- 
vective heat transfer. Only one cover of the beam 
was subjected to heating. It was heated with 12 
quartz-tube, tungsten-filament heating elements uni- 
formly spaced in a reflector. Total output of the 
heaters was 5.7 Btu/see. The reflector, a rectangular 
silver-plated channel the 
heating elements, confronted the heated cover of the 
beam. The reflectance of the reflector surface was 
about 0.97 

Temperatures were measured at the center cross 
section the heated cover of the beam, at the 
center cross section on the unheated cover, at two 
points on the heated cover 1 in. from one end, and at 
one point on the heated cover 1 in. from the other 
end. The temperatures were measured with iron- 
conustantan thermocouples whose outputs were indi- 
eated by galvanometers. The thermocouples were 
calibrated before and after the test. The temper- 
ature rise of the heated cover was approximately 12 
percent the beam ends than at the center 
The experimentally determined heated 
cover temperatures, which were used as the boundary 


Wiis 


stainless-steel enclosing 


On 


less at 


cross section 


condition in eases 1, 2, and 3, are shown in figure 2 
and, together with observed temperatures for two 
points on the opposite cover, in figure 7. 

The normal total emittance of the tvpe 302 
stainless steel was determined by comparing its 


rate of radiation at a given temperature with that of 
a silicon carbide bar at the same temperature. The 
radiation determined by focusing the 


rates were 


images of equal areas of the silicon carbide bar and of 


the stainless steel successively on a thermopile, using 
a fluorite lens. The emittance of the stainless steel 
obtained was constant in the temperature range 400 
to 1,000° F and equal to 0.35. 

Beam deflection at the center cross section 
measured by means of SR-4 type AB—5 electrical 
strain gages mounted back-to-back on a= shielded 
cantilever bye um which Wis deflected by displace- 
ment of the center of the beam relative to its ends 
Contact of the center of the beam with the cantilever 
Wiis made bry il Vveor rod Calibration Wis 
complished by deflecting the cantilever by known 
amounts at the point of contact with a micrometer 


Was 


ac- 


SCTOCW 

The maximum deflection observed was 0.415 in 
after 240 see of heating After 378 see 
the center deflection was 0 37 in. 


of heating 


6. Results 


6.1. Beam-Temperature Distribution 


Computed temperature histories at point «4, in 
the center of the heated cover and at point B, in the 
center of the opposite cover (see sketch on table 1) 
are shown in figures 2 to 6 for the 13 cases considered 
Figure 7 shows the computed temperature distribu- 
tion in a thick-walled beam (cases 1, 2, 3) of the 
same dimensions as the test beam after 148, 189, and 
297 see of heating. Figure 8 shows the temperature 
distribution in a thin-walled beam after 45 see 
slow heating, and after 13.7 see of fast heating. 
Temperatures at A and B are listed in columns 9 
and 10 of table 1. Comparison was made of geo- 
metrically similar beams heated at the same rate 
for the same length of time but having the following 
different values of interior or exterior emittances: 0, 
0.35, and 1.0. 

It was found that, in all cases considered, a change 
in emittance of the surfaces had an appreciable effect 
the temperatures and little on the web 
temperatures when the maximum heated cover 
temperature was above 700° to 900° F. For beams 
with the heated cover at a temperature of about 
1,200° F, changing the internal emittance from 0 to 
1.0 reduced the temperature difference, d, between 
points A and B by 18 percent for the thin-walled 
beam with the cover heated at 91 and by 48 
percent for the same beam with one cover heated at 
30 For the same cover-heating rate the 
effeet of change in interior emittance was less fot 
the thicker-walled beam. 


of 


on cover 


KF sec, 


IK sec. 
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A change in exterior surface emittance has less 
effect on d than ¢ 
tance. Comparison of cases 4 and 6 with corres- 
ponding cases 7 and 8 indicates that, for the particu- 


interior surface emittance from 0 to 1.0 ecauses a 


1 change in interior surface emit- | decrease in d of 19 percent; increasing the exterior 


surface emittance from 0 to 1.0 causes a decrease in 
d of only 8S percent 
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6.2. Thermal Stress and Deflection 


In order to evaluate the effect of a change in 
emittance on the thermal stress distribution and 
deflection, thermal stress distributions and deflec- 
tions were computed for the 13 cases under the 
conditions listed in table 1. Maximum tensile 
stress, Maximum compressive stress, and maximum 
deflection are listed in columns 11, 12, and 13, 
respectively, in the table. To generalize the results 
as regards beam length the deflection is given as the 
deflection per unit beam length squared; the maxi- 
mum deflection ts proportional to the beam length 
squared, the thermal stresses are independent of the 
beam length. The possible effects of vielding, creep, 
and buckling were neglected. Values of the thermal 
stress and deflection are therefore valid only for 
comparison purposes. 

For a thick-walled, slowly heated beam with 
heated cover at 958° F (eases 1, 2, 3) a change in the 
interior wall emittance has little effect on the magni- 
tude of the maximum thermal stress. An increase 
in the interior wall emittance from 0 to 1.0 however, 
results in a decrease of 34 percent in the maximum 
beam deflection. For the beams with the same 
geometry and under the same heating conditions in 
which the final heated cover temperature is about 

200° F, a change of the interior wall emittance 
from 0 to 1.0 reduces the absolute maximum stress 
by from 9 percent (cases 4 and 6) to 27 percent (cases 
9 and 11), and deflection by from 15 percent (cases 
12 and 13) to 40 percent (cases 9 and 11). In the 
latter beam, the reduction in maximum tensile stress 
was 31 percent. For a change of the exterior wall 
emittance from 0 to 1.0 (cases 7 and 8), the redue- 
tion in absolute maximum stress was 8 percent, and 
the reduction in maximum deflection was 7 percent. 

For the beams analvzed, changing the interior 
beam emittance from 0.35 to 1 causes a change in 
the temperature difference, d, between A and B 
which is about 75 percent of the change that would 
occur for a change in emittance from 0 to 1. Cor- 
responding ratios for maximum thermal stress and 
deflection are about 66 and 74 percent, respectively. 





6.3. Comparison of Experimental and 
Theoretical Results 


A partial check of the temperature distribution 
and deflection was made for case 2 using a beam 30 
in. long. Computed and observed temperatures for 
two points on the unheated cover after 148, 189, 
and 297 secs of heating are shown in figure 7. Agree- 
ment is fair. Some of the discrepancy may be 
attributed to gas heat transfer within the beam. 
For the conditions listed under case 2 in table i. 
observed beam deflection was 0.37 in. and beam 
deflection computed from the theoretical temper- 
ature distribution was 0.38 in. i.e., 0.00042 in./(in. 
length)?. 

7. Conclusions 


A change in the interior surface emittance of a 
box beam heated along one cover has an appreciable 
effect on the temperature distribution, and con- 
sequently on the maximum thermal stress and the 
deflection. For a variety of beam geometries and 
heating conditions, a change of the interior surface 
emittance from 0.35 to 1 for a stainless steel beam 
reduces considerably the temperature difference 
from cover to cover in the temperature range above 
700° to 900° F and the maximum thermal stress 
and the deflection in the temperature range above 
900° to 1,200° F. The effect is particularly marked 
for a thin-walled slowly-heated beam. 
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IMPORTANT NOTICE 


Beginning in July 1959, the Journal of Research will be published 
in four separate sections. See the announcement given on the last 
page of this issue. 
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